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Preface 


This  work  extends  the  knowledge  base  for  nonlinear  effects  on  the 
stability  and  control  of  a  spinning  syaaetric  satellite  in  an  elliptical 
orbit.  The  full  nonlinear  controlled  and  uncontrolled  equations  of  aotion 
were  evaluated  using  various  nonlinear  analysis  techniques.  The  origin  of 
an  attractor  in  the  global  phase  space  found  in  Cole’s  work  analyzing  the 
■odal  coordinate  Floquet  control  design  was  the  incentive  for  this 
investigation.  However,  in  order  to  attespt  to  trace  the  origin  of  this 
coaplex  phase  trajectory  structure  in  a  global  sense,  a  siaplified  rate 
and  position  feedback  controller  was  created  and  evaluated.  The  results 
of  this  thesis  show  that  a  global  stability  analysis  yields  results 
particular  to  a  very  specific  set  of  conditions.  One  can  only  hope  to 
characterize  a  ssall  set  of  these  conditions.  But,  there  is  still  nuch  to 
be  learned  about  the  linear  and  nonlinear  stability  characteristics  for 
this  probles. 

I  have  to  thank  Dr.  Wiesel  for  his  advice  and  expertise  on  the 
computer  progress  I  used  for  sany  parts  of  this  study.  Also,  I  must  give 
special  thanks  to  Dr.  Calico  for  the  precious  tine  he  took  fros  his  busy 
schedule  to  assist  se  with  Cole’s  results.  Most  of  all,  I  as  indebted  to 
sy  advisor.  Captain  Jis  Planeaux,  for  his  patience,  insight,  and  guidance 
over  the  past  year.  Finally,  I  as  especially  grateful  to  sy  wife  Patricia 
for  her  love  and  support  that  sade  this  work  possible. 

Michael  Ulisse 
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Abstract 

Controlled  and  uncontrolled  equations  of  aotion  for  a  spinning 
syMetric  satellite  in  an  elliptical  orMt  are  analyzed  using  Floquet 
theory,  Poincare  naps,  and  continuation.  The  equations  of  motion  for  the 
uncontrolled  system  are  derived  via  energy  methods.  A  simple  feedback 
controller  design  and  one  based  on  Floquet  theory  are  presented. 
Parameters  varied  in  the  equations  of  motion  include  the  spin  about  the 
symmetry  axis  of  the  satellite,  the  inertia  ratio  of  the  satellite,  the 
eccentricity  of  the  orbit,  and  the  controller  gains.  Floquet  analyses 
determine  local  stability  and  bifurcation  points.  Poincare  maps  exhibit 
some  interesting  and  complex  structure  in  the  global  phase  space. 
Continuation  using  a  bifurcation  analysis  computer  code  (AUTO)  reveals  a 
sequence  of  period  doublings,  which  is  verified  via  numerical  integration, 
for  the  simple  rate  feedback  controller.  The  onset  of  cheiotic  motion  is 
detected  for  this  sequence  of  period  doublings. 
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NONLINEAR  ANALYSIS  OF  A  SPINNING 


SYMMETRIC  SATELLITE  IN  AN  ELLIPTICAL  ORBIT 


I.  Introduction 

The  attitude  stability  of  a  satellite  under  less  than  noainal 
conditions  is  the  desire  of  any  spacecraft  designer.  Usually,  the  noainal 
conditions  include  operation  of  the  satellite  in  the  linear  regiae  since 
Bost  design  and  stability  analysis  is  easily  acconplished  there.  Yet, 
there  is  always  a  set  of  initial  conditions  that  define  the  limit  of 
effectiveness  of  any  attitude  stabilization  nethod.  These  initial 
conditions  usually  bring  nonlinear  effects  into  the  problea.  Thus,  the 
goal  of  the  satellite  designer  is  to  improve  the  stability  characteristics 
such  that  even  the  less  than  nominal  initial  conditions  pose  no  threat  to 
operational  effectiveness.  Before  such  a  goal  can  be  met,  the  satellite 
stability  characteristics  must  be  traced  from  the  most  basic  design  to  the 
most  complex  to  include  euialyses  away  from  the  linear  regime.  In  other 
words,  a  nonlinear  analysis  may  reveal  some  critical  stability 
characteristics  of  the  system. 

A  passive  method  of  gaining  attitude  stability  is  to  spin  the 
satellite  about  one  of  the  principal  axes  similar  to  a  gyroscope  or  top. 
Kane  nd  Barba  [8]  studied  such  an  approach  for  a  symmetric  satellite  in 
an  elliptical  orbit.  They  found  the  linearized  equations  of  motion  for 
the  satellite  to  be  time  periodic;  Floquet  theory  was  used  to  test  for 
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stability.  They  denonstrated  that  the  attitude  stability  for  a  satellite 
spinning  about  its  syaaetry  axis  is  a  function  of  orbit  eccentricity,  spin 
rate,  and  the  inertia  ratio  of  the  satellite.  Wallace  and  Meirovitch  [16] 
then  investigated  the  stability  characteristics  of  this  problea  using  the 
Liapounov  direct  Method  along  with  asyaptotic  expansions  as  their 
approach.  They  found  that  the  nonlinear  effects  liaited  the  aaplitude  of 
resonant  oscillations  in  the  attitude  Motion  of  the  satellite.  The 
nonlinear  equations  of  Motion  for  this  problea  are  derived  in  Appendix  A 
and  the  linear  equations  are  suMMarized  in  Chapter  I. 

But,  for  fine  pointing  requireaents,  active  control  is  required. 
Siaple  rate  and  position  feedback  aay  be  all  that  is  needed.  However, 
More  coaplex  control  designs  are  possible  and  can  iaprove  perfornance. 
Calico  and  Yeakel  [4]  and  Calico  and  Wlesel  [3]  developed  a  control  Method 
baaed  on  Floquet  theory  utilizing  nodal  coordinates.  In  sinple  terns, 
this  approach  provides  a  tine  periodic,  coupled  rate  and  position  feedback 
to  stabilize  unstable  nodes  (characteristic  Multipliers  outside  the  unit 
circle).  Myers  [12]  evaluated  sone  variations  of  this  approach  while 
Shell  [14]  extended  this  work  to  exanine  what  happens  away  fron  the  linear 
regine  (equilibriun  point).  Cole  [5]  investigated  in  greater  detail  the 
global  effectiveness  of  the  Floquet  controller  by  using  Poincare  naps.  He 
found  an  attractor  present  away  froa  the  noninal  linear  conditions  as  the 
gain  on  the  controller  was  increased.  The  above  work  has  only  begun  to 
scratch  the  surface  of  the  nonlinear  phenonena  present  in  this  problen. 
The  characterization  of  such  phenonena  can  help  define  the  boundary  for 
the  region  of  effectiveness  of  the  Floquet  and  other  controller  designs. 
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The  origin  and  characterization  of  such  coaplex,  global  nonlinear 
phenonena  is  the  goal  of  this  thesis.  To’’ard  this  end,  Floquet  theory, 
Poincare  naps,  bifurcation  theory,  and  continuation  are  used  to  analyze 
controlled  and  uncontrolled  equations  of  notion.  The  uncontrolled 
equations  of  notion  are  analyzed  to  establish  conservative  bifurcation 
points.  The  conplexity  of  the  Floquet  controlled  equations  of  notion  did 
not  allow  for  any  continuation  analyses  on  this  design.  Thus,  this 
controller  design  is  evaluated  in  a  cursory  fashion  in  this  study.  In  its 
place,  a  sinple  rate  and  position  feedback  controller  is  used  to  better 
understeind  local  and  global  stability  characteristics  for  the  controlled 
systen.  The  results  of  this  thesis  add  significantly  to  the  existing 
knowledge  base  for  the  stability  and  control  of  the  syuetric  spinning 
satellite  in  an  elliptical  orbit. 
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II.  Theory 


This  chapter  contains  a  suMary  of  the  equations  of  aotion  for  a 
spinning  syuetric  satellite  in  an  elliptical  orbit,  Floquet  theory,  and 
the  controlled  equations  of  Motion.  The  equations  of  Motion  are  derived 
in  full  detail  in  Appendix  A  via  the  energy  Method  approach.  When  the 
attitude  equations  of  Motion  are  linearized,  they  becoMe  a  set  of  linear, 
tiMe-periodic  ordinary  differential  equations.  This  type  of  systeM  of 
equations  is  essentially  the  Floquet  probleM.  An  extension  of  Floquet 
theory  to  incorporate  Modal  variables  allows  for  the  design  of  a  scalar 
feedback  controller.  A  siMplified  controller  is  then  defined  by  using 
standard  rate  and  position  feedback.  Again,  Floquet  theory  allows  for  the 
deterMination  of  linear  stability  of  both  the  controlled  and  uncontrolled 
equations  of  Motion. 

2.1  Uncontrolled  Equations  of  Motion 

The  equations  of  Motion  for  any  dynanical  systeM  can  be  derived  by 
Newtonian  dynaaics  or  energy  Methods  using  Lagrange’s  equations.  The 
foraer  is  thought  of  as  the  brute  force  approach  for  dynaaical  systeas 
described  by  aany  rotating  coordinate  fraaes  referenced  to  soae  inertial 
fraae.  On  the  other  hand,  the  latter  approach  is  elegant  in  that  the 
equations  of  notion  directly  result  froa  the  expression  for  the  sun  of  the 
kinetic  and  potential  energies  for  the  dynaaical  systen.  In  this  case, 
all  conservative  forces  are  included  as  potentials.  All  nonconservative 
forces  Must  be  included  as  generalized  forces  in  Lagrange’s  equations  of 
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notion.  The  najor  difference  between  the  two  approaches  is  that  the 
energy  nethods  do  not  require  the  inclusion  of  constraint  forces  that  do 
no  work  in  the  problem. 

I  chose  to  present  the  energy  method  approach  in  this  study.  Kane 
and  Barba  [8]  presented  the  Newtonian  derivation  in  their  paper  on  the 
stability  of  the  linearized  attitude  equations  of  notion.  Captain  Dale 
Shell  [14:2.23-2.25]  and  Captain  Jim  Cole  [5:A.1-A.16]  presented  the 
Newtonian  derivations  in  their  theses  also.  A  detailed  presentation  of 
the  energy  nethods  approach  is  provided  in  Appendix  A  for  the  reader.  A 
summary  of  this  is  presented  here. 

The  equations  of  notion  are  for  a  spinning  symmetric  satellite  in  an 
elliptical  orbit  about  a  symmetric  attracting  body.  The  center  of  mass  of 
the  attracting  body  is  considered  fixed.  An  inertial  frame  is  fixed  to 
this  point  and  has  the  I|  and  axes  in  the  equatorial  plane  of  the  body 
with  the  axis  normal  to  this  plane  through  the  north  pole.  The 
satellite  center  of  mass  has  a  coordinate  reference  frame  fixed  to  it  such 
that  axis  Aj  is  in  the  plane  of  the  orbit  and  points  in  the  direction  of 
the  radius  vector  to  the  satellite  center  of  mass.  The  A^  axis  is 
perpendicular  to  axis  A|  and  in  the  direction  of  the  satellite’s  orbital 
notion.  Finally,  the  Aj  axis  is  simply  in  the  direction  which  completes 
a  dextral,  orthogonal  coordinate  system.  This  puts  the  Aj  cucis  normal  to 
the  orbit  plane.  A  1-2-3  rotation  sequence  through  the  angles  6p  62,  and 
Bj  measure  the  attitude  orientation  of  the  satellite  with  respect  to  the 
orbital  reference  frame  A.  Using  these  definitions,  the  equations  of 
notion  in  nondimens ional  coordinates  are 
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C"  C’ 


(2.1) 


o'/- 2®i0a  t  autjBj + 2  v'BaCOsdj +¥''00861  tan6j  +  v'^cosOja  indj 
-(jr+  lXs+l)(6a  •*■  v'sin6i)sec6j 


(2.2) 


0^'-(#r+ lXfl+ 1)(6iC086j  -  v'cosOiSinOj) 

-|6l*-v'*C08*6i+^j8ln62C0862-v''8in6i-2v'6'iC086iC0a*6, 


(2.3) 


where 


v-iizsifii 

c* 


(2.4) 


C’ 


(2.5) 


6j  -(8+  l)t  -V  -(8+  1)J1 1-¥ 


(2.6) 


'’'he  diaensionless  orbital  radius  is  defined  as 


(2.7) 


with  nondiaensional  tiae  aiaply  as 


T»J3t 


(2.8) 


which  is  the  aean  action  of  the  satellite  tiaes  tiae.  The  inertia  ratio 
and  the  spin  parameter  are  defined  by 


^9 


(2.9) 
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respectively.  Note  that  the  orbital  aotion  equation  given  by  equation 
(2.1)  drives  the  attitude  equations  to  be  tiae  periodic.  See  Appendix  A 
for  definitions  of  all  other  variables. 

By  inspection,  a  constant  solution  to  the  attitude  equations  exists 
and  is  given  by 

di-da-O  (2.11) 

This  equilibrium  solution  keeps  the  spin  axis  of  the  satellite 
perpendicular  to  the  orbit  plane.  If  the  attitude  equations  are 
linearized  about  this  equilibrium  point,  we  get  a  set  of  linear,  time- 
periodic,  ordinary  differential  equations.  In  order  to  study  the 
stability  of  these  attitude  equations,  we  must  use  Floqu6t  theory. 

2. 2  FloquSt  Theory 

Picking  a  state  vector  as 

jp-{eie,e'ei}  (2.12) 

we  can  write  the  system  of  nonlinear  equation  as 

(2.13) 

The  system  of  equations  can  be  linearized  by 

(2.14) 


where 


(2.15) 


d£ 

W 


:A(t) 


Note  that  the  A  aatrix  is  evaluated  at  the  equilibriua  condition.  The  A 
■atrix  can  be  written  as 


At)= 


*11 

*12 

*13 

*14 

*21 

*22 

*23 

*24 

*31 

*32 

*33 

*34 

*41 

*42 

*43 

*44 

(2.16) 


The  resulting  teras  of  the  A  aatrix  for  the  linearized  systea  are 


(2.17) 


(2.18) 


*34 


(2.19) 


*42-*3l'^ 


(2.20) 


a^3-«a4-l  (2.21) 

with  the  reaaining  teraa  being  equal  to  zero.  The  orbital  aotion  of  the 
satellite  causes 
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(2.22) 


In  other  words,  the  A  aatrix  is  periodic  with  a  period  of  T.  We  now  have 
a  set  of  linear,  tiae  periodic,  ordinary  differential  equations.  Floquet 
was  the  first  to  describe  the  solution  to  such  a  problea  during  the  last 
decade  of  the  19th  century  [3:671], 

The  solution  to  the  this  type  of  system  of  equations  can  be  shown  to 
be 

(2.23) 

where  4(t,0)  is  the  state  transition  aatrix.  The  state  transition  aatrix 
obeys  the  relationships 

(2.24) 

and 

♦(0,0)-J  (2.25) 

where  I  is  the  identity  matrix  [3:671].  Now,  from  Floquet  theory,  the 
state  transition  matrix  becomes 

(2.26) 

where  J  is  the  constant  diagonal  matrix  of  Poincar^/characteristic 
exponents  auid  F(t)  is  a  periodic  matrix  with  same  period  T  as  the  original 
system  of  equations.  Therefore,  the  periodicity  of  the  F  matrix  allows  us 
to  say  that 

F{0)mf{T)  (2.27) 

at  time  T.  This  reduces  equation  (2.26)  to 
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^T.0)^F{0)e'"F-\0)  (2.28) 

4(T,0)  is  an  iaportant  Batrix;!!  is  called  the  Bonodrosy  Batrix. 
NuBerical  integration  over  the  period  T  of  equation  (2.24)  with  initial 
conditions  given  by  equation  (2.25)  will  deteraine  the  BonodroBy  Batrix. 

Since  the  F  Batrix  is  periodic,  it  is  bounded  and  does  not  affect 
the  stability  of  the  systeB.  So,  the  only  way  for  the  systeB  to  be 
unstable  is  for  the  unbounded  growth  of  the  exponential  terB.  Since  J  is 
a  diagonal  Batrix,  it  follows  that  the  F  Batrix  is  the  Batrix  of 
eigenvectors  for  the  BonodroBy  Batrix.  So,  the  eigenvalues  of  the 
BonodroBy  Batrix  are  related  to  the  Poincare/characteristic  exponents  by 

(2.29) 

where  the  ^^’s  are  the  Poincare/characteriatic  exponents  and  the  A-’s  are 
the  eigenvalues  of  the  BonodroBy  Batrix.  The  eigenvalues  are  called  the 
Floqu4t/characteristic  aultipliers.  So,  if  the  Boduli  of  all  the 
characteristic  aultipliers  are  less  than  one,  the  systea  is  linearly 
stable  about  the  equilibriua  condition.  If  any  of  the  eigenvalues  have 
Boduli  greater  than  one,  then  the  systea  is  unstable.  If  any  of  the 
eigenvalues  have  Boduli  equal  to  one,  then  the  systea  is  only  stable  if 
the  Bultiplicity  of  that  particular  root  is  equal  to  the  nullity  of  the 
BonodroBy  aatrix  [8:404]. 

2.3  Controlled  Equations  of  Motion 

The  controlled  attitude  equations  of  aotion  are 
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Oi-2O'i02tan62+2v^92CO86i+v''cosditan02+v'’coB6jSin0i  (2  30) 
-(JC+  IXs*  1)(0^ + v'sin6i)sec6j + jq0i+ +  Jl^0i + 

03 1XS+ 1)(0iCO802  -  v'coB0i  B  in6j) 

-|^,’-v»cos^,»^jslnO,co.ej  12.31) 

-v"Bin0i-2v'0ico8OiCO8^j 

+Ai0i+^0j+^^i+lt;0'2 


The  additional  ter«8  added  to  equations  (2.30)  and  (2.31)  provide  the 
necessary  feedback  to  cancel  disturbances  away  fros  the  desired 
equilibriuB  condition  for  an  otherwise  unstable  uncontrolled  systea. 
Floqu6t  theory  and  the  use  of  sodal  coordinates  allow  for  the  design  of 
this  tise  periodic  controller.  This  particular  design  approach  aoves  one 
pair  of  unstable  characteristic  aultipliers  into  the  unit  circle.  In 
other  words,  the  controller  forced  the  characteristic  aultipliers  to  a 
stable  regiae  of  the  phase  space.  This  technique  was  developed  by  Calico 
and  Wiesel  [3:673-674].  The  reader  is  also  referred  to  [5]  and  [14]  for 
further  reading  on  this  developaent. 

The  in  equations  (2.30)  and  (2.31)  are 


X^*-OiXi3i+<2jI.j3 


(2.32) 


where  and  Gj  are  gains  that  deteraine  the  location  of  the  characteristic 
aultipliers  in  the  coaplex  plane.  The  are  found  by  solving  the  linear, 
tine-periodic  aatrix  differential  equation 
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L'(t)-I<t)4j-A  *([t)X<t)  (2.33) 

where  A(t)  is  given  by  equation  (2.16)  and  J  is  the  constant  aatrix  of 
characteristic  exponents.  This  equation  is  a  reoult  of  aanipulating 


equation  (2.26)  froB  Floquet  theory  where 

L^{t)^F-\t)  (2.34) 

Again,  see  [3]  and  [5]  for  further  detail. 

A  simplified  controller  is  defined  by  the  following  equations: 

Oi-26^02tan62+2v'0aco86i+v''co86itan6j+v'’cos6i8in6ji  (2  35) 
-<jr+lX«+l)(6a +v'8in6j8ec62+Gidi+G^0j  *0^9(1*0^91^ 

9(3  "(K*  IX^*  1)(6iC0862' Vco80^8in62) 

~y"e  InBj  -  2  v'8ico80iCO8*6j 

♦  Gj®  ♦  Gj^B2  ♦  ©40^  ♦  0^03 

Now,  the  G^’s  in  equation  (2.36)  are  simple  feedback  gains  which  are 
arranged  such  that  there  can  be  uncoupled  position  and  rate  feedback. 
These  equations  are  easier  to  work  with  since  the  complexities  of  the 
Ploqu6t  controller  are  not  included. 

Floquet  theory  is  used  to  determine  the  stability  of  both  controlled 
systems  of  equations.  Numerical  integration  and  Poincar6  plots  are  used 
to  study  the  stability  of  the  local  and  global  solutions.  One  goal  will 
be  to  use  these  techniques  to  facilitate  the  design  of  a  control  for  this 
problem. 
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III.  Numerical  Methods  for  Nonlinear  Analysis 


This  chapter  provides  a  brief  sumaary  of  a  few  selected  nuaerical 
nonlinear  analysis  techniques.  The  discussion  starts  with  an  introduction 
to  the  nonlinear  ordinary  differential  equations  problem.  The  computer’s 
ability  to  help  investigate  such  problems  is  highlighted.  The  thrust  of 
the  remaining  nonlinear  analysis  discussion  emphasizes  bifurcation  theory. 
Bifurcations  due  to  characteristic  multiplier  location  and  movement  and 
their  effect  on  stability  are  described.  The  next  two  sections  summarize 
the  selected  numerical  techniques  for  nonlinear  analysis.  These 
techniques  include  Poincare  maps  and  continuation  using  AUTO  [6]  which  is 
a  computer  code  used  for  part  of  this  study.  For  completeness,  a 
discussion  of  Hamiltonian  (conservative)  systems  is  included.  Specific 
problems  encountered  with  AUTO  for  Hamiltonian  systems  are  included  in  the 
discussion. 

3.1  Nonlinear  Analysis  and  the  Computer 

In  many  fields  of  science  and  engineering,  dynamics  of  motion  or 
processes  are  modelled  by  ordinary  differential  equations.  In  most  cases, 
these  models  contain  nonlinearities  such  as  equations  (2.2)  and  (2.3).  If 
there  is  no  explicit  dependence  on  time  in  the  differential  equation,  then 
the  system  is  said  to  be  autonomous.  A  general,  second  order,  autonomous 
differential  equation  can  be  written  as 
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where  the  vector  x  contains  the  states  of  the  system  and  f  is  the  vector 
of  functions  describing  the  time  evolution  of  the  vector  x. 

Since  equation  (3.1)  is  a  nonlinear  set  of  differential  equations, 
one  would  expect  nonlinear  phenomena  in  the  solutions.  However,  the 
solutions  to  such  problems  are  rarely  obtained  in  closed  form. 
Traditionally,  the  equations  are  linearized  about  some  equilibrium  point 
which  gives  a  local  picture  of  the  what’s  going  on  in  the  dynamics  of  a 
particular  system.  The  objective  of  such  an  approach  is  to  stay  away  fro* 
the  nonlinearities  of  the  problem  at  hand.  In  reality,  we  can  not  be 
assured  of  initial  conditions  that  keep  us  in  the  local  (linear)  regime. 
Thus,  a  global  picture  away  from  the  equilibrium  point  is  required  in 
order  to  address  the  nature  of  any  nonlinear  phenomena  in  the  problem. 

With  the  advent  of  the  computer  and  efficient  and  accurate  numerical 
integration  methods,  the  study  of  nonlinear  phenomena  has  become  possible 
for  those  with  little  expertise  in  more  rigorous,  mathematical  approaches. 
Using  this  straightforward  approach,  one  can  begin  to  see  the 
nonlinearities  effect  on  the  time  evolution  of  a  dynamical  system. 
"Jumping  between  modes,  sudden  onset  or  vanishing  of  periodic  solutions, 
loss  or  gain  of  stability,  buckling  of  frames  and  shells,  ignition  and 
combustion,  auid  chaos  are  but  a  few  examples”  [13:  ix].  So,  through  the 
use  of  the  computer,  numerical  techniques  have  opened  the  door  to  the 
study  of  nonlinear  phenomena. 


3.2  Bifurcation  of  Solutions 


The  introduction  of  a  free  parameter  into  the  differential  equation 
modelling  the  dynamics  of  system  can  cause  a  variety  of  responses  for 
certain  values  of  the  parameter.  The  dynamical  model  is  now  represented 
by 

where  1  is  a  free  parameter.  A  bifurcation  occurs  when  a  change  in 
response  type  is  detected  as  the  parameter  passes  through  a  critical  value 
called  the  bifurcation  point  [7:317].  A  bifurcation  diagram,  which  plots 
the  value  of  the  varying  parameter  versus  some  scalar  measure  of  one  of 
the  system  states,  allows  us  to  study  the  effect  of  the  parameter  on  the 
system  response.  In  other  words,  bifurcation  of  the  system  response  can 
be  detected  by  looking  at  the  bifurcation  diagram.  Much  of  this  section 
follows  the  presentation  in  Seydel  [13]  which  is  an  excellent  reference 
for  bifurcation  theory. 

The  characteristic  multipliers  of  a  periodic  system  determine  its 
local  (linear)  stability.  Recall  from  the  discussion  on  Floquet  theory, 
that  the  moduli  of  the  characteristic  multipliers  directly  indicate 
whether  or  not  there  is  local  stability.  For  a  periodic  system,  there 
will  be  one  characteristic  multiplier  with  a  modulus  of  unity.  If  all  of 
the  other  characteristic  multipliers  lie  in  the  unit  circle,  the  periodic 
motion  is  stable.  In  order  for  local  stability  to  be  lost,  one  or  more  of 
the  characteristic  multipliers  must  depart  the  confines  of  the  unit 
circle.  This  movement  of  the  characteristic  multipliers  is  a  result  of 
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varying  the  free  parameter  X.  Note  that  soae  dynaaical  systeas  aay  have 
a  vector  of  paraaeters  that  can  be  varied.  Figure  3.1,  borrowed  froa 
oeydel  [13:299],  depicts  three  possible  ways  of  leaving  the  confines  of 
the  unit  circle. 


In  Figure  2.1  (A),  the  characteristic  aultiplier  crosses  the  unit 
circle  at  tl.  This  type  of  crossing  results  in  one  of  three 
possibilities.  One  possibility  is  that  a  liait  point  is  achieved.  In 
other  words,  for  solutions  exist  with  one  being  stable  and 
the  other  being  unstable.  If  ,  then  no  solutions  exist.  The  other 
two  possibilities  are  a  result  of  soae  fora  of  syaaetry  in  the  problea 
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which  results  in  pitchfork  or  transcritical  bifurcations  [13:2521. 

In  Figure  3.1  (B),  the  characteristic  aultiplier  crosses  the  unit 
circle  at  -1.  This  type  of  crossing  usually  results  in  a  period  doubling 
for  the  dynaaical  systes.  This  type  of  bifurcation  is  only  possible  for 
systeas  with  three  dimensions  or  higher.  This  type  of  bifurcation, 
soaetiaes  called  a  flip  bifurc&tion,  exhibits  a  loss  of  stability  for  the 
original  stable  periodic  orbit  while  a  new  orbit  with  twice  the  period  is 
created  [15:127]. 

In  Figure  3.1  (C),  a  pair  of  complex-conjugate  multipliers  cross  the 
unit  circle  for  a  critical  value  of  the  parameter  1.  This  type  of 
bifurcation  is  known  as  a  bifurcation  into  a  torus.  Again,  such  a 
bifurcation  only  exits  for  three  dimensions  or  higher.  The  original 
periodic  orbit  loses  stability  and  remains  at  the  center  of  the  torus. 
These  tori  can  be  attractive  or  repelling  [13:263].  Figure  3.2  depicts 
motion  on  a  torus.  A  slice  of  the  torus  taken  with  a  Poincare  map  yields 
a  closed  curve  on  the  nap. 

Motion  on  a  torus  contains  two  frequencies.  One  frequency  is  for 
notion  around  the  circumference  of  the  torus  and  the  other  frequency  is 
for  notion  around  the  cross-section  of  the  torus.  If  these  two 
frequencies  are  irrational,  then  the  notion  is  said  to  be  quasiperiodic. 

Two  less  understood  but  most  interesting  cases  of  bifurcation  are 
those  that  result  in  a  strcuige  attractor  or  chaos.  An  attractor  is  called 
str  ..ige  if  the  trajectory  depends  sensitively  on  initial  conditions 
[13:310].  A  definition  for  a  strange  attractor  is  "an  attractor  that  is 
not  an  equilibrium  point  nor  a  limit  cycle,  nor  a  quasiperiodic  attractor" 
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Figure  3.2  Motion  on  a  Torus 

[11:267].  Chaos  is  not  far  reaoved  froB  a  strange  attractor.  Both  are 
sensitive  to  initial  conditions.  However,  trajectories  starting  fros 
slightly  different  initial  conditions  exponentially  diverge  for  chaotic 
■otion  [11:262].  It  is  alsost  absurd  to  think  that  a  detersinistic 
differential  equation  results  in  an  unpredictable  type  of  probles  (chaos). 

3,3  PoincarS  Maps 

A  technique  which  is  attributed  to  Henri  Poincare  (1854-1912)  is  the 
Poincar4  sap.  It  is  also  known  in  other  terns  as  a  first  return  sap, 
surface  of  section,  and  Poincar4  section.  "Poincare  sapping  is  often  used 
for  analyzing  the  geometry  of  trajectories  in  the  phase  apace"  [9:26]. 
For  trajectories  in  the  phase  8p{u:e  that  are  periodic,  the  Poincar6  sap  is 
the  logical  choice  to  sake  to  identify  the  structure  of  any  attractors. 


3-6 


An  attractor  is  a  portion  of  phase  space  that  trajectories  have  an 
affinity  to  after  the  transients  have  all  died  out.  Classical  dynaaical 
attractors  include  equilibriun  points,  lisit  cycles,  and  toroidal  surfaces 
[11:261]. 

The  Poincare  sap  is  the  result  of  a  trajectory  in  the  phase  space  as 
it  repeatedly  passes  transversely  through  a  plane  or  surface  that  slices 
through  the  phase  space.  The  sequence  of  points  plotted  on  the  surface 
forms  the  nap.  The  characteristic  structure  of  the  Poincar4  nap  can  be 
used  to  describe  the  attractor  at  hand.  Table  3.1  derived  directly  froa 
Moon  [11:53]  sunsarizes  vrhat  particular  structures  on  the  Poincare  nap 
night  indicate. 


Table  3.1  Classification  of  Poincare  Maps 


finif  number  ot  pointai  periodic  or  subhanMnic  occlllatlon 
Cloaad  curv«:  quaalpar Iodic,  two  locnnwmraf  frosuMclM  protMiit 
Open  curva:  auggamt  •odolllng  as  a  on#  dlaanalonal  aap;  try  plotting  x(t)  varaaa  x(t«T) 
fractal  collaction  of  polnta:  atranga  attractor  In  tbraa  pltaaa  apaca  diaanalona 
fuzzy  collaction  ot  pointai  (1)  dynaaical  ayataa  with  too  auch  randoa  or  nolay  Input; 
(11)  strange  attractor  but  ayataa  haa  vary  aaall  dlaslpatlon;  (111)  atranga  attractor 
In  phaaa  space  with  aora  than  three  diaanalona;  (Iv)  quaalparlodlc  action  with  threw  or 
aora  doalnant  Incoaaanaurata  fraquanclaa 

[ll:53J 


A  variation  of  the  Poincar4  nap  is  the  stroboscopic  technique 
[15:66].  This  approach  creates  a  Poincar4  nap  by  plotting  points  of  the 
trajectory  in  the  phase  space  at  none  specified  tine  which  is  usually  a 
nultiple  of  the  period  of  the  systen. 

In  higher  dinensional  systens,  it  is  difficult  to  identify  the 
characteristics  of  the  attractor  in  the  phase  space  using  the  standard 
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Poincar^  aap.  The  problea  here  is  trying  to  view  aore  than  three 
disensions  on  a  two  diaensional  plot.  A  technique  called  the  double 
Poincare  map  allows  one  to  observe  the  structure  of  a  three  diaensional 
Poincare  sap  of  a  fourth  order  systea  [11:142].  The  double  Poincare  map 
is  constructed  by  saapling  the  phase  space  trajectory  as  done  with  the 
stroboscopic  technique  ;however,  two  frequencies  are  used  instead  of  one. 
A  point  is  plotted  "when  the  first  frequency  is  coincident  with  the  second 
frequency"  [11:143].  This  technique  cuts  a  slice  of  the  three  diaensional 
phase  space  in  order  to  deteraine  the  nature  of  any  attractors  present. 
Any  fractal  nature  in  the  resulting  plot  is  characteristic  of  a  strange 
attractor  [11:144]. 

3.4  Hamiltonian  Systems 

A  Haailtonian  systea  is  synonymous  with  a  conservative  systea;  no 
dissipation  is  present.  Many  classical  djniaaical  probleas  such  as  those 
encountered  in  orbital/celestial  aechanics  fall  into  this  category  of 
probleas.  These  types  of  probleas  are  characterized  by  forces  that  can  be 
represented  by  a  potential.  The  uncontrolled  equations  of  aotion  in  this 
thesis  are  conservative  since  the  only  force  in  the  problea,  gravity,  is 
represented  by  potential  as  shown  in  Appendix  A. 

A  Haailtonian  systea  will  always  have  an  even  nuaber  of  phase 
variables.  If  a  Poincare  plot  is  used  to  study  the  trajectories  of  a  four 
state  Haailtonian  systea,  the  phase  trajectories  will  eventually  pass 
through  all  points  in  the  reaaining  three-diaensional  subspace.  "Chaotic 
orbits  in  conservative  systeas  tend  to  visit  all  parts  of  a  subspace 
uniforaly;  that  is  they  exhibit  a  unifora  probability  density  over 
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restricted  regions  in  the  phase  space"  [11:63].  It  is  important  to  note 
that  conservative  systess  "do  not  exhibit  strange  attractors"  [13:311]. 

If  dissipation  is  added  to  a  conservative  system  by  some  form  of 
control,  one  would  expect  a  contraction  of  the  trajectories  in  the  phase 
space  into  some  form  of  defined  structure.  Chaotic  trajectories  in 
dissipative  systems  exhibit  a  fractal  structure  in  the  phase  space 
[11:63].  However,  a  small  amount  of  dissipation  in  a  system  can  cause  it 
to  resemble  a  Hamiltonian  system  in  the  phase  space  [11:63].  So,  it 
follows  that  the  addition  of  enough  dissipation  should  result  in  the 
formation  of  some  sort  of  attractor. 

3.5  Continuation 

A  numerical  method  that  allows  us  to  compute  a  family  of  solutions 
(algebraic  or  periodic)  as  a  parameter  is  varied  given  an  initial  solution 
is  a  continuation  method.  Another  definition  is  "a  numerical  technique 
enabling  us  to  obtain  one  branch  of  solutions  (or  more  branches  of 
solutions  mutually  connected  at  branch  points)"  [9:21].  Continuation  is 
also  known  as  branch  tracing.  Sc,  with  continuation,  a  parameter  is 
systematically  varied  to  compute  new  solutions  from  nearby  previously 
computed  solutions.  Stability  of  these  solutions  are  available  since  the 
Jacobian  matrix  for  the  system  is  also  carried  along  in  the  solution 
process.  Therefore,  it  is  simply  a  matter  of  finding  the  eigenvalues  of 
the  Jacobian  matrix  to  determine  stability.  As  the  eigenvalues  cross  the 
unit  circle,  bifurcation  points  are  located. 

A  FORTRAN  digital  computer  program  created  by  E.  Doedel  called  AUTO 
[6]  performs  continuation  analysis  for  moderate  sized  problems.  It 
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coaputes  faailies  of  solutions  for  both  equilibria  and  periodic  solutions. 
This  software  package  was  used  in  this  study  in  a  liaited  fashion  due  to 
the  conservative  nature  of  probleas  studied  herein.  AUTO  was  designed  to 
study  dissipative  systeas.  In  order  to  handle  conservative  (Haailtonian) 
types  of  probleas,  dissipation  aust  be  added  with  caution  so  that  the 
problea  being  studied  is  not  altered  greatly. 

Take  for  exaaple,  equation  (2.1),  which  is  the  nondiaensional 
differential  equation  for  the  satellite’s  orbital  aotion.  Since  this 
equation  is  a  result  of  a  conservative  systea,  an  infinite  faaily  of 
closed  trajectories  exist.  Such  a  systea  causes  a  problea  in  AUTO.  So, 
dissipation  aust  be  added  in  order  for  AUTO  to  run  correctly  and  converge. 
In  other  words,  a  liait  cycle  about  the  desired  orbit  aust  be  created  by 
the  addition  of  a  dissipative  tera  to  equation  (2.1).  The  question  that 
reaains  is  how  to  accoaplish  such  a  feat. 

A  chapter  in  [1]  is  devoted  to  such  a  topic.  Let  us  rewrite 
equation  (2.1)  as  set  of  first  order  differential  equations. 

(3.3) 


(3.4) 


H  can  be  thought  of  as  the  Haailtonian  for  the  systea.  The  eliaination  of 
the  nondiaensional  tine  variable  gives  us  the  equation  for  the  phase  paths 
of  the  systea  which  is 


dy.  if  (!-•*)  1 

dx  y[  jc*. 


(3.5) 
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Since  equation  (3.5)  is  separable,  it  is  integrated  to  get 


=  (3.6) 

2  2x^  X 

where  H  is  the  constant  which  describes  a  particular  orbit.  For  the 
nondimensional  variables  used,  H.  is  equal  to  for  all  possible  orbits 
with  initial  conditions  at  apogee  of  y=0  and  x=l+e.  However,  we  are 
interested  in  the  particular  case  e=0.5.  The  feature  we  want  the 
dissipation  term  to  have  is  a  zero  contribution  when  the  trajectory  is  on 
the  desired  orbit  and  an  other  than  zero  contribution  elsewhere.  A 
dissipative  term  of  the  form 

tf{x,y)[H^-Hix,y)]  (3.7) 

will  cleverly  do  this  for  us.  Now,  if  we  let  the  new  system  of  equations 
with  the  added  dissipation  be 

,3.8, 


(3.9) 

w!u-re  5  is  a  parameter,  we  can  show  that  a  limit  cycle  about  the  desired 
orbit  is  created  by  tlie  use  of  the  invariant  integral 

j" j  t:i[J]cixdy=0  (3.10) 

'  t-  -b:.'.:  iUi  ;.iaTr:\  ;’c  r  t;;e  system.  It  can  ije  si;o;.n  ‘hat  the 
iii’-arruit  inti.-gral  holds  for  equations  (3.8)  and  (3.9)  over  the  region  of 
the  desired  orbit.  Also,  it  can  be  shown  that  the  limit  cycle  created 
about  the  desired  orbit  is  stable. 
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IV.  Results 


This  chapter  contains  a  variety  of  results  obtained  fro*  using  the 
theory  and  techniques  described  in  chapters  II  and  III.  The  objective  of 
such  a  broad  approach  was  to  gain  enough  information  about  the  syuetric 
satellite  in  an  elliptical  orbit  problem  to  characterize  some  of  the  local 
and  global  phenomena  as  a  function  of  parameters  in  the  system.  The 
parameters  varied  included  the  spin  rate  about  the  axis  of  symmetry  (s), 
inertia  ratio  (K),  eccentricity  (e),  and  controller  gains  (Gj,  Gj,  G^). 
Both  global  and  local  stability  characteristics  were  examined.  However, 
the  local  stability  characterization  was  a  more  tractable  problem.  As  you 
will  see,  the  nonlinearities  in  this  problem  created  some  interesting 
results  which  were  at  times  difficult  to  explain. 

The  computer  was  the  primary  tool  used  in  generating  the  results 
presented  in  this  thesis.  All  the  computer  code  used  was  written  in 
FORTRAN  77  using  double  precision  variables.  The  routines  used  for  the 
Floquet  amalysis  originated  from  previous  work  [14],  [12],  and  [5].  I 
modified  these  routines  to  allow  for  the  cases  studied  in  this  thesis  as 
well  as  automating  the  program  to  vary  a  chosen  parameter  over  a  specified 
range  at  a  user  selected  step  size.  These  routines  were  run  on  a  Digital 
Equipment  Corporation  (DEC)  VAX  8550  computer  using  the  VAX  operating 
system.  The  numerical  integration  of  the  equations  of  motion  was 
performed  using  the  Hamming  integrator.  This  integration  was  done  on  the 
Elxsi  6400  computer  using  the  EMBOS  operating  system.  The  AUTO 
continuation  routines  were  run  on  the  DEC  VAX  8550. 
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4.1  Floquet  Analysis 

Floquet  theory  was  used  to  analyze  characteristic  Multiplier 
Movement  on  the  complex  plane  to  deteraine  possible  bifurcation  points. 
It  was  hoped  that  any  bifurcation  points  found  would  lead  to  soae  of  the 
complex  structure  found  by  Cole  [5].  Local  stability  characterization  in 
a  finer  detail  than  presented  by  Kane  and  Barba  [8]  resulted  from  this 
analysis. 

4.1.1  Uncontrolled  Equations  of  Motion.  Kane  and  Barba  [8],  Myers 

[12],  and  Shell  [14]  all  calculated  characteristic  multipliers  for  the 

uncontrolled  equations  of  motion.  However,  these  multipliers  were  only 

calculated  at  values  of  eccentricity  equal  to  0.0,  0.1,  0.3,  0.5  and  0.7 

for  a  set  of  chosen  values  of  the  spin  parameter  and  the  inertia  ratio. 

In  this  thesis,  the  characteristic  multipliers  for  this  problem  were 

-8 

calculated  over  the  range  O.OleSO.?  with  steps  in  e  as  small  as  1  x  10 
to  determine  bifurcation  points.  The  bifurcation  points  for  the 
uncontrolled  equations  are  termed  conservative  because  the  multipliers 
split  away  from  the  unit  circle  in  pairs.  These  multipliers  were  plotted 
on  the  complex  plane  along  with  the  unit  circle  to  give  a  finer  picture  of 
what  is  happening  as  e  is  varied. 

The  first  case  examined  contained  the  parameter  values  used  in  the 
Floqu4t  controller  design  and  analyses  done  in  the  past  [12],  [14],  [5]. 
Figure  4.1  shows  this  most  interesting  case  with  8=1.0  and  K=0.7.  At 
e=0.5,  which  is  the  case  the  Floqu4t  controller  stabilized,  one  multiplier 
pair  is  outside  the  unit  circle  and  one  multiplier  pair  is  inside  the  unit 
circle.  Since  one  pair  of  multipliers  is  outside  the  unit  circle,  there 
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is  no  local  stability  for  the  equilibriun  condition.  Figure  4.1  exhibits 
three  places  where  the  aultipliers  leave  the  unit  circle.  One  place  is 
represented  by  a  dark  filled  in  oval  on  the  unit  circle.  Figures  4.2, 
4.3,  4.4,  and  4.5  detail  the  values  of  a  where  the  Bultipliers  leave  the 
unit  circle;  these  points  locate  where  local  stability  is  lost.  Figure 

4.2  shows  a  plot  of  the  Multipliers  for  0.38<e<0.41.  This  departure  from 
the  circle  was  completely  unexpected  and  not  found  previously.  Figure 

4.3  shows  that  the  multipliers  leave  the  unit  circle  for  0.46<e<0.47. 
When  a  complex  pair  of  multipliers  leave  the  unit  circle  as  in  figures  4.2 
and  4.3,  it  usually  indicates  a  bifurcation  into  a  torus.  Figure  4.4 
shows  that  the  multipliers  return  to  the  unit  circle  for  0. 6316<eS0. 6317 
where  local  stability  is  regained  for  a  very  small  range  of  e.  This 
region  of  local  stability  was  not  found  previously.  Figure  4.5  shows  how 
sensitive  the  multipliers  can  be  to  a  very  small  change  in  e.  For 
0.63274217<e<0. 63474218,  the  multipliers  leave  the  unit  circle  with  a 
bifurcation  point  that  resembles  a  local  period  doubling. 

The  effect  of  the  spin  parameter  on  figure  4.1  was  looked  at  next. 
Figures  4.6,  4.7,  4.8,  4.9,  4.10,  4.11,  4.12,  and  4.13  show  what  happens 
to  local  stability  at  different  values  of  the  spin  parameter  s  in  the 
range  -3.0<8<3.0. 

First,  the  spin  parameter  was  increased  over  the  value  used  in 
figure  4.1.  Figures  4.6  and  4.7  show  the  movement  of  the  characteristic 
multipliers  for  the  spin  parameters  8=1.5  and  8=2.0  re8i>ectively.  The 
detailed  structure  of  leaving  and  returning  to  the  unit  circle  has 
disappeared.  However,  the  period  doubling  bifurcation  point  still 
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Uncontrolled  System  Multipliers 
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Uncontrolled  System  Multipliers 
(0.632742  1  7<e<0.632742 1 8.s=  1 .0,K=0  7) 


reaains.  This  occurs  between  0. 547<e<0. 548  in  Figure  4.6  and  between 
0. 506<e<0. 507  in  Figure  4.7.  An  explanation  for  the  greater  stability  in 
this  case  is  that  the  increased  spin  has  added  to  the  angular  aoBentua, 
thus  increasing  stability.  Figures  4.8  and  4.9  show  a  further  increase  of 
the  spin  parameter  to  the  values  of  3=2.5  and  8=3.0  respectively.  There 
is  still  a  period  doubling  bifurcation  for  8=2.5  between  0.699<e^.7. 
However,  for  s=3.0,  there  is  no  loss  of  stability;  all  characteristic 
multipliers  remained  on  the  unit  circle  for  0.0<e<0.7.  Again,  this  is 
expected  for  high  spin  rates  about  the  83rmaetry  axis  of  the  satellite. 

Next,  the  spin  parameter  was  reduced  from  the  baiseline  value  of 
8=1.0  used  in  figure  4.1.  For  s=0.0,  the  rotation  of  the  satellite  about 
the  Dj  axis  is  simply  equal  to  the  mean  motion  of  the  satellite.  From 
Figure  4.10,  we  see  two  bifurcation  points  present.  One  is  a  period 
doubling  type  occurring  between  0.25<ei0.251  amd  the  other  is  a  pitch  fork 
bifurcation  or  limit  point  between  0.490<e<0.491.  One  pair  of  multipliers 
return  back  to  the  unit  circle  for  0.601Se<0.602.  Figure  4.11  shows  the 
characteristic  multipliers  for  s=-1.0.  One  multiplier  is  off  the  graph 
making  the  system  locally  unstable  for  0.0<e<.0.7.  This  agrees  with  Sana 
and  Barba  [8:404-405]  and  makes  sense  since  the  net  spin  of  the  satellite 
is  zero  here.  Figure  4.12  is  imother  interesting  case  with  characteristic 
multipliers  leaving  and  returning  to  the  unit  circle  with  s=-2.0.  The 
multipliers  leave  and  return  to  the  unit  circle  for  0.496<e^. 533.  Again, 
this  is  a  local  bifurcation  into  a  torus.  There  is  a  pitchfork  type 
bifurcatior.  or  limit  point  for  0.677Se<0.678  also.  This  structure 
disappears  for  s=-3.0.  However,  figure  4.13  shows  a  period  doubling 
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bifurcations  for  0. 49<e<0. 491 .  Again,  with  the  larger  aagnitude  spin 
parameter  in  these  two  cases,  local  stability  is  maintained  over  a  larger 
region. 

The  effect  of  the  inertia  ratio  was  studied  next.  Figures  4.14, 
4.15,  4.16,  4.17,  4.18,  4.19,  4.20,  and  4.21  show  the  effect  of  the 
inertia  parameter  K  for  0.0<e<0.7.  The  spin  parameter  was  kept  at  3=1.0 
like  the  baseline  case  in  figure  4.1  A  positive  K  indicates  an  oblate 
satellite  configuration  while  a  negative  K  indicates  a  prolate 
configuration. 

First,  the  oblate  values  for  K  were  examined.  Figure  4.14  shows  the 
characteristic  multipliers  for  K=1.0  for  0.0<e<0.7.  There  is  a  period 
doubling  loss  of  stability  for  0.268<e<0.27  and  0.58^^.582.  Again,  the 
complex  structure  of  the  baseline  case  in  figure  4.1  is  not  present 
anymore.  Figure  4.15  shows  the  characteristic  multiplier  movement  for 
K=0.5  in  which  the  multipliers  all  stay  on  the  unit  circle  for  0.0<e<0.7. 
This  agrees  with  Kane  and  Barba  except  for  e=0.7  where  they  gave  a  loss  of 
local  stability  at  K=0.4  and  K=0.6;  one  would  infer  that  K=0,5  is  also 
unstable.  Figure  4.16  shows  a  bifurcation  into  a  torus  for  0. 552<eS0. 554 
with  K=0.25.  Figure  4.17  indicates  that  an  inertia  parameter  of  0.0 
remains  locally  stable  for  0.0le<0.7.  For  such  a  case,  the  gravitational 
potential  given  by  equation  (A. 58)  reduces  simply  to  the  potential  acting 
on  the  center  of  mass  of  the  system.  In  other  words,  there  are  no  torques 
on  the  satellite  due  to  gravity. 

The  prolate  satellite  configurations  were  studied  next.  Figures 
4.18  and  4.19  exhibit  some  very  interesting  movement  of  the  characteristic 
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■ultipliers  at  K=-0.25  and  K=-0.5  respectively.  At  K=-0.25,  two 
■ultipliers  aeet  at  +1  on  the  unit  circle  between  0. 362Se<0. 364  putting 
all  the  Bultipliers  on  the  unit  circle.  Two  coaplex  pairs  of  multipliers 
bifurcate  into  a  torus  for  0. 434<e<0.436  with  one  set  of  multipliers  going 
inside  the  unit  circle  and  one  set  going  outside.  At  K=-0.5,  we  note  a 
period  doubling  bifurcation  between  0.240^^.242  and  a  pitchfork  tyi)e 
bifurcation  or  limit  point  between  0.05e<0.002.  Note  how  the  multipliers 
move  around  the  outside  of  the  unit  circle  in  figure  4.19.  Figures  4.20 
and  4.21  have  only  one  pair  of  characteristic  multiplier  inside  the  unit 
circle  for  K=-0.75  and  K=-1.0  respectively.  The  other  two  multipliers  in 
these  cases  are  outside  the  unit  circle  and  do  not  appear  on  these  graphs. 
This  indicates  local  instability  for  0.0Se<0.7  which  agrees  with  Kane  and 
Barba. 

In  summary,  a  more  complex  and  detailed  structure  of  local  stability 
characteristics  using  plots  of  characteristic  multipliers  for  the  various 
conditions  considered  herein  was  revealed.  This  is  certainly  more 
information  than  previously  appreciated  for  the  spinning  symmetric 
satellite  in  an  elliptical  orbit  problem. 

4.1,2  Floquet  Controlled  Equations  of  Motion.  The  Ploquet 
controller  was  designed  to  cancel  disturbances  away  from  the  desired 
equilibrium  condition  with  the  values  of  the  parameters  being  esO.5, 
8=1.0,  K=0.7.  Figure  4.22  shows  how  one  pair  of  characteristic 
multipliers  is  moved  inside  the  unit  circle  as  the  Floqu4t  controller 
feedback  gains  G1  and  G2  in  equations  (2.30)  and  (2.31)  are  varied  between 
0.0  and  1.0.  The  other  pair  is  not  affected  since  this  is  a  "modal 
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control"  design  approach.  There  is  only  a  torus  bifurcation  point  on  this 
plot  of  characteristic  aultipliers.  Cole  claiaed  that  a  period  doubling 
bifurcation  existed  in  a  plot  he  presented  of  the  characteristic  exponents 
[5:4-15-4-16].  However,  figure  4.22  does  not  indicate  such  a  period 
doubling  bifurcation.  Even  though  Cole  was  in  error  here,  he  still  say 
have  found  a  global  period  doubling  as  his  results  sees  to  suggest.  The 
only  conclusion  that  can  be  sade  froa  figure  4.22  is  that  local  stability 
is  achieved  by  the  noveaent  of  the  unstable  pair  of  aultipliers  into  the 
unit  circle. 

4.1.3  Ra.te  and  Position  Feedback  Equations  of  Motion.  Rather  than 
varying  the  paraaeters  in  the  Floquet  controlled  equations  of  action,  a 
siaple  rate  and  position  feedback  controller  was  designed.  This  was  done 
to  allow  for  the  use  of  the  continuation  software  (AUTO)  and  is  explained 
later.  In  order  to  save  tiae  in  intialising  the  continuation  software,  I 
perforaed  a  Floquet  analysis  on  the  siaple  controller  design  to  locate  any 
existing  bifurcation  points.  Once  any  bifurcation  points  were  found,  it 
was  easy  to  start  the  continuation  software  up  with  the  correct  initial 
conditions  for  paraaeters  in  the  equations  of  action. 

First,  only  uncoupled  rate  feedback  was  added  to  the  equations  of 
action  given  by  equations  (2.35)  and  (2.36).  The  gain  Gj  was  varied  for 
0.0<Gj:S1.0  to  achieve  local  stability.  Figure  4.23  indicates  how  effective 
this  siaple  controller  was  in  aoving  the  unstable  pair  of  aultipliers  into 
the  unit  circle.  In  coaparison  with  figure  4.22,  we  see  a  siailar 
aoveaent  of  the  unstable  pair  of  aultipliers.  But,  as  exhibited  by  figure 
4.23,  the  siaple  controller  also  aoved  the  stable  pair  of  aultipliers 
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further  into  the  unit  circle.  The  Floquet  controller  left  this 
Bultiplier’s  location  unchanged.  Note  that  there  is  a  torus  bifurcation 
point  as  a  coaplex  pair  of  nultipliers  cross  the  unit  circle. 

The  uncoupled  rate  feedback  gain  was  fixed  to  Gj=0.85  so  that  all  the 
characteristic  aultipliers  would  just  lie  on  the  real  axis  in  figure  4.23. 
This  condition  of  local  stability  was  now  analyzed  as  other  parameters  in 
the  equations  were  varied.  Figure  4.24  plots  the  BoveBent  of  the 
characteristic  Bultipliers  as  the  inertia  ratio  went  froB  -1.0<K<1.0. 
When  K  was  varied  through  this  range,  two  crossings  of  the  unit  circle  at 
+1  occurred.  These  pitchfork  type  bifurcations  happened  in  the  ranges  of 
-0. 36<K<-0. 35  and  -0. 29SK<-0. 28.  Both  of  these  crossings  Boved  a 
Bultiplier  into  the  unit  circle  for  a  positive  change  in  K.  The  other  two 
Bultipliers  stayed  in  the  unit  circle  over  the  range  K  was  varied.  So, 
after  the  second  crossing  into  the  unit  circle,  local  stability  was 
achieved.  We  can  conclude  that  the  simple  rate  feedback  controller  with 
G^=0.85  is  locally  stable  for  -0.28<K<1.0.  The  second  crossing  was 
analyzed  using  the  continuation  software  with  the  results  presented  later 
in  this  thesis. 

Figure  4.25  depicts  characteristic  multiplier  movement  as  the  spin 
parameter  is  varied  from  -3.0<823.0.  As  with  the  previous  case,  two 
multipliers  remain  inside  the  unit  circle  and  two  start  outside  the  unit 
circle.  Again,  there  are  two  crossings  into  the  unit  circle  at  +1.  One 
occurred  for  -1.82isi-1.81  and  the  other  between  -0.23SaS-0.22.  These  two 
bifurcations  are  pitchfork  type  bifurcations.  All  the  multipliers  entered 
inside  the  unit  circle  after  the  second  crossing,  thus  indicating  local 
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stability  for  -0.22<sS3.0. 

Next  the  eccentricity  was  varied  for  0.0<e<0.7.  Figure  4.26  shows 
that  for  s=1.0,  K=0.7  and  Gj=0.85,  local  stability  is  maintained  over  the 
total  range  of  eccentricity.  So,  the  simple  rate  feedback  eliminated  any 
effect  that  eccentricity  might  have  on  this  particular  case. 

Finally,  coupled  rate  feedback  and  position  feedback  were  added  to 
the  equations  of  motion  given  by  equations  (2.35)  and  (2.36).  These  gains 
are  and  Gj  respectively  in  the  equations  of  motion.  They  were  both 
varied  individually  between  0.0  and  1.0.  Again,  the  fixed  parameters  were 
Gj=0.85,  e=0.5,  s=1.0,  K=0.7.  Figures  4.27  and  4.28  both  indicate  that 
coupled  rate  feedback  by  itself  and  coupled  position  feedback  by  itself 
have  little  effect  on  the  stability  of  the  system. 

Uncoupled  rate  feedback  with  GjsO.85  was  chosen  to  b<;  the  baseline 
design  for  fixed  parameters  of  e=0.5,  s=1.0,  KaO.7.  Figures  4.29,  4.30, 
4.31  ,4.32  show  that  this  gain  selection  returns  the  system  to  equilibrium 
from  an  initial  disturbance  of  0j=02=O.4.  Note  that  it  took  between  two 
and  three  orbits  to  cancel  the  disturbance. 

4,2  Poincare  Maps 

Poincare  maps  were  plotted  for  the  uncontrolled,  Ploquet 
controlled,  and  rate  feedback  controlled  equations  of  motion.  These  maps 
were  plotted  to  give  a  snapshot  of  the  global  phase  space  at  a  variety  of 
initial  values  for  the  parameters  in  the  equations  of  motion.  It  was 
hoped  that  some  of  the  cases  selected  would  shed  some  light  on  the  origin 
of  the  attractor  found  by  Cole  in  [5].  The  stroboscopic  technique  for 
creating  Poincare  maps  was  used.  Most  cases  plotted  a  point  on  the  map  at 
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a  2n  interval  since  the  orbital  period  is  2n.  But,  soae  cases  used  a  4ii 
interval  instead.  The  technique  produced  saps  depicting  soae  interesting 
and  coaplex  structure  in  the  gJobal  phase  space. 

4.2.1  Uncontrolled  Equations  of  Motion.  Two  cases  were  looked  at 
for  the  uncontrolled  equations  of  aotion.  One  case  was  locally  stable 
with  all  of  the  characteristic  aultipliers  on  the  unit  circle  with  e=0.15. 
The  other  case  was  locally  unstable  with  e=0.5.  These  two  cases  were 
chosen  for  two  basic  reasons.  One  is  that  the  stable  case  occurs  before 
the  torus  bifurcation  point  and  the  unstable  case  occurs  after  the  torus 
bifurcation  point.  This  bifurcation  point  is  shown  by  figure  4.3.  The 
second  reason  is  that  the  Floqu4t  controller  was  studied  at  e^O.S.  The 
hypothesis  here  is  that  the  global  phase  space  structure  found  by  Cole 
originated  fros  the  uncontrolled  systes.  The  saps  plotted  2500  orbits  at 
a  2ji  interval  with  fixed  parameter  values  of  s=1.0  and  K=0.7.  Both  cases 
used  a  disturbance  of  6js62=0.005  to  create  the  aaps. 

The  Poincar^  naps  for  the  es0.15  case  exhibit  an  identifiable 
structure.  Figure  4.33  shows  a  aap  for  0j  vs  @2  and  figure  4.34  shows  a 
nap  for  the  rates.  The  trajectories  seen  to  have  two  frequencies  plus  the 
orbital  frequency  which  is  characteristic  of  notion  on  a  3-torus.  This 
aakes  sense  since  only  a  snail  step  was  taken  away  fron  the  locally  stable 
periodic  orbit  and  both  sets  of  aultipliers  are  coaplex  with  a  nodulus  of 
unity.  The  structure  in  both  of  these  aaps  is  confined  to  a  region  in  the 
phase  space. 

On  the  other  hand,  the  e=0.5  case  shows  no  identifiable  structure  on 
the  Poincar^  aaps.  These  aaps  are  shown  in  figures  4.35  and  4.36.  Note 
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that  the  naps  are  still  confined  to  a  region  of  the  phase  space.  This 
result  is  typical  of  a  Haailtonian  systea  since  the  aap  fills  the  space 
occupied  without  any  preferred  structure.  If  daaping  is  added  to  the 
systea,  one  expects  this  region  to  collapse  into  an  attractor  though. 

4,2.2  Floquet  Controlled  Equations  of  Motion.  A  few  selected 
Poincar4  aaps  were  generated  for  the  Floquet  controlled  equations  of 
motion.  The  paraaeter  values,  which  were  identical  to  the  values  used  by 
Cole,  were  selected  as  e=0.5,  8=1.0,  K=0.7,  and  Gj=G2=0.4.  These  aaps  were 
plotted  to  see  if  any  additional  inforaation  could  be  gained.  The  initial 
disturbance  put  into  the  systea  was  0j=O.O851,  0j=O.O576,  0j=-O.1569,  and 

t 

02=0.95.  In  aodal  variables  used  by  Cole,  this  was  equivalent  to 
q|=i|2=q]=0.0  and  q^=0.95.  All  the  aaps  in  this  case  were  saapled  at  a  An 
interval  since  Cole  claiaed  the  systea  had  period  doubled  [5]. 

The  attractor  found  by  Cole  is  recreated  here.  Figures  4.37  and 
4.38  show  the  attractor  which  was  created  by  the  first  400  orbits  aapped 
in  the  aodal  variables.  These  aaps  contain  two  apparently  closed  curves 
instead  of  four  as  deaonstrated  by  Cole  in  [5].  The  shapes  in  the  figures 
natch  Cole’s  identically.  The  reason  for  Cole  having  four  closed  curves 
is  that  he  actually  plotted  at  a  2k  interval  rather  than  4k. 

When  2500  orbits  are  napped,  a  different  picture  is  created. 
Figures  4.39  and  4.40  show  the  original  closed  curves  eabedded  in  sone 
other  structure  in  the  phase  space.  There  are  now  also  what  seens  to  be 
solid  lines  appearing  on  the  aaps.  The  points  which  created  these  solid 
lines  caae  froa  the  very  last  portion  of  the  2500  points  plotted.  It 
appears  that  the  initial  closed  curves  are  transients  that  disappear  over 
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the  long  time  scale  and  maybe  re-visited  over  an  even  longer  time  scale. 

The  same  type  of  structure  appeared  in  the  nondimensional  variable 
phase  space.  Figures  4.41  and  4.42  map  the  nondimensional  physical 
variables  using  the  same  parameter  values  as  for  the  modal  variables. 
However,  these  figures  map  an  angle  versus  its  corresponding  rate  for  2500 
orbits.  Again,  there  are  two  closed  curves  embedded  in  some  other 
structure  in  the  phase  space.  Also,  the  solid  looking  lines  are  present. 

The  solid  lines  seemed  to  indicate  that  another  attractor  was 
present.  It  is  possible  that  this  attractor  is  a  limit  cycle  that  is 
wandering  a  little  due  to  numerical  accuracy  problems  in  the  integration. 
So,  the  integration  was  restarted  using  the  end  points  of  figures  4.41  and 
4.42  for  2500  more  orbits.  The  results  of  this  integration  are  mapped  in 
figures  4.43  and  4.44.  Figure  4.43  shows  two  different  closed  curves 
embedded  in  some  other  structure  in  the  phase  space.  Again,  the  solid 
lines  appear.  Figure  4.44  is  a  little  harder  to  visualize.  However,  two 
closed  curves  are  present  inside  the  structure  on  the  map.  These  curves 
occurred  during  the  first  part  of  the  mapping.  The  same  solid  lines 
appear  again.  If  the  first  2500  orbit  and  the  2500  restarted  orbits  are 
combined  onto  one  map,  we  see  four  closed  curves.  This  is  easily  seen  in 
figure  4.45.  Yet,  there  is  still  other  structure  present  on  the  map. 
Notice  that  the  solid  lines  fall  on  top  of  one  another.  The  same 
phenomenon  occurs  in  figure  4.46.  These  maps  show  that  the  trajectory 
visits  two  weakly  unstable  structures  (one  like  a  torus  and  the  other  like 
a  limit  cycle)  with  transients  between  the  two  leading  to  a  complicated 
attractor.  An  explanation  of  this  may  be  that  the  trajectories  in  the 
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phase  space  are  trapped  between  the  surface  of  an  unstable  torus  and  an 
unstable  lisit  cycle.  The  solid  lines  on  the  plot  may  indicate  the 
approach  of  the  trajectories  to  the  li«it  cycle. 

4.2.3  Rate  and  Position  Feedback  Equations  of  Motion.  A  small 
amount  of  damping  was  added  to  the  uncontrolled  equations  of  motion  here 
through  the  uncoupled  rate  feedback  term.  Again,  the  case  we  are 
interested  in  here  is  when  e=0.5.  Recall  that  the  uncontrolled  equations 
of  motion  showed  no  structure  on  the  Poincare  map  at  this  particular 
eccentricity.  See  figures  4.35  and  4.36.  With  the  added  damping,  we  were 
looking  for  the  formation  of  an  attractor  similar  in  structure  to  that  of 
the  Floquet  controlled  equations.  But,  there  was  no  guarantee  that  this 
would  be  the  case. 

The  problem  here  was  set  up  by  plotting  the  characteristic 
multipliers  for  the  rtte  feedback  gains  of  GjsQ.O,  Gj=0.01,  and  G2=0.03  for 
0.46<e<0.5,  s=l.€),  K=0.7.  Figure  4.47  shows  the  plot  for  these  three  rate 
feedback  gains.  The  solid  line  in  the  plot  is  a  portion  of  the  unit 
circle.  The  location  of  the  characteristic  multipliers  for  e=0.5  is  the 
first  point  on  left  of  each  plotted  line.  For  G^=0.01,  the  multipliers  are 
still  outside  the  unit  circle.  It  can  be  deduced  from  figure  4.47  that 
all  the  multipliers  are  inside  the  unit  circle  for  Gj=0.05  and  e=0.5. 
These  two  cases  were  looked  at  with  the  Poincare  map. 

The  maps  for  Gj=0.01  indicate  an  initial  transient  attractor  as  shown 
earlier.  Figures  4.48,  4.49,  4.50,  and  4.51  show  this  for  500  orbits. 
The  disturbance  for  figure  4.48  and  4.49  was  6j=02=O.OO5.  The  initial 
disturbance  for  figures  4.50  and  4.51  was  0|=92=O.4.  Note  that  the  same 
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attractor  appeared  for  two  disturbances  differing  by  two  orders  of 
■agnitude  and  two  sampling  intervals  of  2ii  and  4n  respectively. 

When  10000  orbits  were  mapped  using  the  smaller  of  the  two 
disturbances,  an  interesting  structure  appeared  in  the  phase  space. 
Figures  4.52  and  4.53  depict  the  original  attractor  embedded  in  some 
intricate  structure.  An  explanation  for  this  is  that  the  added 
dissipation  was  so  small  that  the  phase  space  still  "wants  to"  resemble 
that  of  the  original  Hamiltonian  system.  But,  some  form  of  structure  is 
present. 

As  Gj  was  increased,  the  intricate  structure  present  began  to  shrink 
in  size  on  the  maps.  When  02=0.05,  an  attr8M:tor  appeared  at  the  origin 
which  was  without  structure.  This  limit  cycle  attractor  is  exhibited  by 
figures  4.54  and  4.55  by  the  fixed  point  at  the  origin  of  each  map.  This 
result  was  what  we  expected  the  rate  feedback  to  accomplish.  The  larger 
of  the  two  disturbances  was  used  here.  So,  the  rate  feedback  term  did  not 
preserve  to  the  global  attractor  present  for  the  Floqu^t  controlled 
system.  This  is  a  more  attractive  approach  to  obtain  the  desired  control 
since  the  global  structure  is  not  present  for  larger  disturbances  away 
from  the  equilibrium  condition. 

4.3  Continuation 

Continuation  was  used  in  a  limited  fashion  in  this  thesis  because  of 
the  conservative  nature  of  the  symmetric  spinning  satellite  problem.  The 
continuation  software,  AUTO,  primarily  handles  dissipative  systems  of 
equations.  This  is  because  conservative  systems  have  multiple  multipliers 
at  unity  which  causes  a  non-uniqueness  of  closed  orbits.  The  software  was 


4-65 


tricked  into  working  for  the  uncontrolled  equations  of  aotion  by  aodifying 
the  orbit  equations  of  notion.  The  software  never  worked  for  the  Floquet 
controlled  equations  of  notion  because  of  the  abundance  of  nultipliers  at 
unity.  However,  the  software  did  work  for  the  rate  and  position  feedback 
equations  of  notion.  For  these  equations,  nore  than  one  branch  of 
solutions  were  calculated  fron  an  initial  bifurcation  point  found  during 
the  Floquet  analyses  presented  earlier. 

4.3.1  Modified  Orbit  Equations  of  Motion  for  AUTO.  The  orbit 
equation  given  by  equation  (2.1)  has  an  infinite  fanily  of  closed-orbit 
solutions.  This  non-uniqueness  was  the  crux  of  the  problen  for  AUTO. 
Figure  4.56  shows  a  fanily  of  solutions  for  various  initial  conditions  on 
a  phase  portrait  for  e=0.5.  A  program  created  by  Dr.  Viiesel  was  used  to 
create  this  and  other  phase  portraits  presented  here. 

The  problen  (i.e.  trick)  here  was  to  create  a  stable  linit  cycle  at 
a  desired  orbit  of  interest.  All  other  closed  curves  in  the  conservative 
phase  portrait  of  (2.1)  must  be  made  to  disappear  to  resolve  the  non¬ 
uniqueness  problen.  One  of  the  desired  orbits  is  shown  in  figure  4.57  for 
e=0.5.  This  orbit  has  a  perigee  of  (=0.5  and  an  apogee  of  (=1.5.  To 
achieve  this  desired  orbit  fron  any  initial  condition,  a  dissipation  tern 
must  be  added  to  the  equation  of  notion  for  the  orbit.  This  was  described 
in  chapter  III,  section  3.5.  The  dissipation  tern  is  shown  in  equation 
(3.8).  Its  nain  feature  is  that  it  is  zero  on  the  orbit  of  interest  and 
not  equal  to  zero  elsewhere. 

A  phase  portrait  of  the  modified  equation  of  notion  for  the  orbit 
exhibited  a  linit  cycle  at  the  chosen  orbit  fron  the  original  conservative 
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system.  Figure  4.58  shows  how  trajectories  starting  inside  and  outside 
the  desired  orbit  converge  on  the  created  limit  cycle.  This  addition  of 
dissipation  allowed  AUTO  to  work  for  the  uncontrolled  and  rate  feedback 
equations. 

4.3.2  Uncontrolled  Equations  of  Motion.  Using  the  modified 
dissipative  orbit  equations,  the  continuation  software  was  run  in  an 
attempt  to  locate  bifurcation  points  as  parameters  were  varied.  Yet, 
since  the  attitude  equations  were  still  conservative,  no  bifurcation 
points  were  detected  because  crossings  occurred  in  multiple  pairs  while 
AUTO  was  looking  for  a  single  pair  to  cross.  However,  the  software  did 
compute  the  characteristic  multipliers  correctly.  Figure  4.59  shows  a 
plot  of  characteristic  multipliers  with  the  parameter  values  used  to 
generate  figure  4.1.  Both  plots  contain  the  same  information.  Note  that 
AUTO  also  included  the  orbit  equation  characteristic  multipliers.  One 
multiplier  is  a  +1  while  the  other  moves  on  the  real  axis  inside  the  unit 
circle.  This  shows  that  the  dissipation  term  has  indeed  created  a  limit 
cycle.  This  method  of  obtaining  the  characteristic  multipliers  was  not 
efficient.  So,  few  runs  were  made  for  this  purpose. 

4.3.3  Floquet  Controlled  Equations  of  Motion.  Nonautonomous 
systems  are  not  well  suited  for  analysis  using  AUTO.  Usually,  auxiliary 
differential  equations  must  be  added  to  the  original  system  of  equations 
to  make  the  total  system  autonomous.  These  auxiliary  differential 
equations  should  have  solutions  that  allow  for  the  calculation  of  the 
nonautonomous  variable.  This  method  of  tricking  the  software  seemed 
attractive  but  was  very  difficult  to  implement  without  causing  ot..  ir 
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problems.  The  original  Floquet  controlled  equations  of  motion  fit  into 
this  category  of  problems  since  the  approximation  of  the  periodic  feedback 
was  given  by  a  time  dependent  Fourier  series.  No  simple  way  was  found  to 
make  the  Floquet  equations  of  motion  autonomous. 

First,  the  linear,  time-periodic  matrix  set  of  differential 
equations  defining  the  periodic  control  was  appended  to  the  system  of 
equations  to  provide  the  controller  feedback,  thus  eliminating  the  need 
for  the  fourier  series.  This  matrix  set  of  equations,  given  by  equation 
(2.33),  added  sixteen  auxiliary  differential  equations  to  the  original  set 
of  six  which  described  the  Floquet  controlled  system.  The  continuation 
software  never  converged  for  this  approach.  This  set  of  equations  was 
then  numerically  integrated  using  Hamming  to  look  for  possible  errors.  It 
was  found  that  the  controller  performed  as  expected  over  a  short  time 
scale.  But,  over  the  long  time  scale,  the  states  of  the  auxiliary 
differential  equations  diverged.  Since  the  controller  equations  were 
derived  basically  from  the  equations  of  variation,  there  was  no  guarantee 
that  they  would  be  stable.  So,  a  Floquet  analysis  was  performed  on  the 
controller  differential  equations.  This  generated  a  16x16  monodromy 
matrix  resulting  in  sixteen  characteristic  multipliers.  Table  4.1  lists 
these  characteristic  multipliers  and  their  moduli.  There  are  four  pairs 
on  the  unit  circle,  two  pairs  outside  the  unit  circle,  and  two  pairs 
inside  the  unit  circle.  This  verifies  that  these  equations  are  locally 
unstable.  It  can  be  concluded  that  the  four  pairs  of  multipliers  at  unity 
will  lead  to  the  non-uniqueness  problem  for  AUTO  discussed  earlier.  No 
other  valid  conclusions  were  reached  on  why  AUTO  did  not  converge. 
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Table  4.1  Characteristic  Multipliers  for  Control  Generation  Equations 


MULTIPLIER 

1 multiplier! 

-0. 129753±1. 242112 

1.248881 

1.248881 

1.248881 

1.248881 

1.248881 

-0.10389510.994588 

1.000000 

-0.10389510.994588 

1.000000 

1.000000 

1.000000 

1 . 000000 

1.000000 

1 . 000000 

1 . 000000 

1.000000 

1.000000 

-0.00831910.796383 

0.800717 

0.800717 

0.800717 

0.800717 

0.800717 

Finally,  sixteen  auxiliary  differential  equations  were  added  to  the 
original  systea  of  differential  equations  whose  solutions  provided  the 
first  eight  trigonoaetric  teras  of  the  fourier  series  used  to  approxiaate 
the  control.  In  other  words,  these  auxiliary  differential  equations  had 
solutions  of  sin(t),  co8(t)  and  the  higher  order  haraonics  which  aade  the 
systea  of  equations  autonoaous.  But,  the  convergence  problea  still 
existed.  It  was  hypothesized  that  the  haraonics  of  the  auxiliary 
equations  caused  the  problea  for  AUTO.  So,  in  the  end,  the  continuation 
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software  could  not  be  used  to  evaluate  the  Floquet  controlled  equations  of 
motion. 

4.3.4  Rate  and  Position  Feedback  Equations  of  Motion.  There  were 
two  plots,  figures  4.24  and  4.25,  of  characteristic  multipliers  for  the 
rate  feedback  controller  that  exhibited  bifurcation  points  that  crossed 
the  unit  circle  at  +1.  The  bifurcation  point  found  in  figure  4.24  was 
analyzed  using  the  continuation  software.  The  results  that  generated 
figure  4.24  located  a  bifurcation  point  between  -0.29<K<-0.28.  AUTO  was 
initialized  to  compute  branch  1  which  is  the  parameter  dependent  branch  of 
trivial  solutions  that  contains  this  bifurcation  point.  AUTO  located  a 
pitchfork  type  bifurcation  point  at  K=-0. 2837217.  The  results  from  AUTO 
agreed  with  the  Floquet  analysis  in  that  the  solution  became  stable  at 
this  bifurcation  point  as  the  multiplier  crossed  into  the  unit  circle. 

AUTO  was  then  used  to  continue  from  this  bifurcation  point  to 
generate  branch  2  of  the  solution.  Since  branch  2  was  continued  from  a 
pitchfork  bifurcation,  two  symmetric  branches  exist  that  are  non-trivial. 
This  branch  exhibited  a  typical  period  doubling  bifurcation.  That  is  that 
the  solution  became  unstable  after  the  bifurcation  at  the  critical 
parameter  value  with  the  birth  of  a  stable  orbit  with  twice  the  period. 
Successive  branches  up  to  branch  5  were  computed  with  the  continuation 
software  to  reveal  a  sequence  of  period  doublings.  Figures  4.60,  4.61, 
4.62,  and  4.63  plot  the  bifurcation  diagrams  using  6^  and  62  maximum  as 
measure  of  the  solution  response  for  these  results.  Table  4.2  then 
summarizes  the  bifurcation  points  found  by  AUTO. 

A  sequence  of  period  doublings,  if  it  continues  ad  infinitum. 
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Bifurcation  Diagrom  of  K  vs  0,  MAX  for  Rote  Feedback  Systeu 
(e^O  ‘^j,s=  I  0,03-0  85) 


150 


Branches  3  and  4  of  Bifurcation  Diagram  of  K  vs  0|  MAX 
fur  Rate  Feedback  System  (e-0  5,s=1  O.G,-0.8S) 


5  and  4  of  Bifurcotion  Diagram  of  K  vs  Qp 
eedbock  System  (e  =  0  5.s- 1 .0,G»=  0  FfS) 


Table  4.2  SuMary  of  Bifurcation  Points 


BRANCH 

TYPE 

Nthitical 

MAXIMUM  6j 

MAXIMUM  02 

1 

Pitchfork 

-0.2837217 

0.000000 

0.000000 

2 

Period  Double 

-0.3358687 

0.967213 

0.871232 

3 

Period  Double 

-0.3438680 

1.061325 

1.051574 

4 

Period  Double 

-0.3440923 

1.075763 

1.060851 

5 

Period  Double 

-0.3451118 

1.077220 

1.065944 

satisfies  the  relationship  discovered  by  Fiegenbaua  which  is 

-4.6692  (4.1) 

The  A’s  in  this  equation  are  the  critical  paraaeter  values  at  which  the 
bifurcation  takes  place.  The  Fiegenbaua  sequence  for  the  critical 
bifurcation  paraaeters  of  this  case  is 

a^-7. 734023, 4. 936993,-  (i-l,»)  (4.2) 

It  appears  that  this  sequence  of  period  doubling  nuabers  is  approaching 
the  Fiegenbaua  constant. 

The  tiae  histories  of  0j  and  02  at  different  values  of  K  depicted 
the  effect  of  the  period  doubling  on  the  attitude  aotion.  Figures  4.64 
and  4.65  are  the  tiae  histories  provided  by  AUTO  for  K=-0. 339954.  This 
point  occurs  on  branch  3  of  the  solutions.  Figures  4.66  and  4.67  are  the 
tiae  histories  obtained  by  nuaerical  integration  between  orbit  499  and 
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-0.50 


-0.50 


501.  The  time  histories  for  both  AUTO  and  numerical  integration  are 
identical  except  for  a  small  phase  shift.  The  time  histories  show  smooth, 
continuous  motion  with  a  4n  period.  A  Poincare  map  was  also  plotted  for 
this  case.  The  map,  shown  in  figure  4. 68,  contains  two  fixed  points  which 
is  expected  for  a  period  doubled  orbit.  So,  the  three  nonlinear  analysis 
tools  used  in  this  thesis  all  reached  the  same  conclusion  here. 

As  an  aside,  the  phase  shift  difference  between  the  data  from  AUTO 
and  numerical  integration  was  looked  into.  A  Poincare  map  for  2500  orbits 
was  generated  with  points  being  plotted  at  4ii  intervals.  One  would  expect 
a  single  point  on  the  map.  However,  figure  4.69  shows  a  line  instead. 
So,  a  time  history  was  plotted  for  data  between  orbits  2001  and  2003  in  an 
attempted  to  find  the  problem.  This  plot  is  shown  in  figure  4.70.  When 
compared  to  figure  4.66,  we  see  a  very  small  phase  shift  between  the  two. 
Evidently,  there  must  be  some  error  in  the  numerical  integration  that 
slowly  causes  the  two  Poincare  map  fixed  points  to  shift  over  the  long 
time  scale.  This  explains  the  line  on  the  Poincare  map. 

In  theory,  a  sequence  of  period  doublings  usually  results  in  onset 
of  chaotic  motion.  So,  the  parameter  value  was  set  to  K=-0.3485  (well 
beyond  the  last  doubling  in  the  sequence)  to  see  if  this  was  the  case 
here.  The  time  histories  between  orbit  200  and  300  are  shown  in  figures 
4.71  and  4.72.  The  motion  now  looks  less  regular  with  no  periodicity 
detectable.  A  Poincare  map  was  plotted  for  this  case  also  to  see  if  any 
identifiable  structure  appears  in  the  phase  space.  Figure  4.73  shows  the 
map  for  1000  orbits  sampled  at  a  2n  interval.  It  does  not  show  any  well 
defined  structure.  So,  it  appears  that  this  sequence  of  period  doublings 
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FiLjure  4  68  F-'oifirare  Mol.  of  Gi  v5:  Gj  from  Numerical  Inti-gration 

(60  Orbits  -lOrnpled  ot  =  0  b,s- 1  0,K=  0,.3  vjyb4,G^= 


Iiift  (e=0  5,5.-  1 .0,l<  -  0  .559954,03-0  85) 


mericnily  Integroted  Orbits  of 


leads  to  chaotic  aotion. 

In  summary,  the  uncoupled  rate  feedback  control  with  a  gain  of 
Gj=0.85  was  relatively  robust  in  that  it  did  not  exhibit  a  global  attractor 
for  larger  disturbances.  This  simple  approach  seemed  to  work  well  for  a 
variety  of  conditions.  However,  extreme  cases  led  to  a  sequence  of  period 
doublings . 
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y.  Conclusions  and  Recommendations 


The  local  stability  analysis  provides  a  finer  detailed  picture  of 
the  linear  stability  characteristics  for  the  uncontrolled  and  controlled 
equations  of  notion  than  was  known  previously.  The  effect  of  the 
paraneters  varied  show  the  nany  possible  solutions  for  the  attitude  notion 
of  the  systems  considered  here.  The  existence  of  conservative  period 
doubling,  pitchfork  or  Unit  point,  and  torus  bifurcations  show  the 
importance  of  analyzing  the  uncontrolled  system  in  detail  first.  The  sane 
type  of  analyses  for  the  controlled  equations  yields  the  sane  type  of 
critical  control  design  infornation.  Basically,  local  stability  is 
characterized  for  nominal  and  less  than  noninal  conditions.  These  results 
establish  a  range  of  operating  conditions  near  the  equilibrium  solution. 
A  rate  feedback  gain  of  0.85  is  reconnended  over  a  wide  range  of 
conditions. 

The  global  stability  analysis  indicates  both  expected  and  unexpected 
structure  in  the  phase  space.  For  example,  the  collapse  of  a  region  of 
phase  space  into  an  attractor  with  addition  of  damping  through  a  range  of 
feedback  gain  is  expected.  But,  a  torus-like  attractor  for  the  Floquet 
controller  is  not  easy  to  explain.  The  nonlinear  analyses  do  bound  the 
Units  of  operation  for  the  systems  studied.  However,  each  case  studied 
is  particular  to  the  parameters  and  the  initial  conditions  selected. 
Broad  statements  of  global  stability  characteristics  should  be  made  with 
extreme  caution.  The  characterization  of  nonlinear  effects  is  indeed  a 
difficult  task  to  accomplish.  Evidently  though,  simple  rate  feedback 
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might  introduce  less  global  behavior  than  Floquet  control. 

The  continuation  analysis  shows  some  of  the  hidden  phenomena  of  this 
system.  The  presence  of  a  period  doubling  may  be  of  interest  if  such 
motion  is  detected  in  an  operational  design.  This  type  of  analyses  gives 
the  control  designer  a  valuable  tool  to  characterize  stability.  However, 
the  continuation  software  used  here  does  not  work  well  for  conservative 
systems  such  as  the  one  studied  here  without  some  careful  modification. 

The  global  stability  analysis  for  this  problem  was  too  time 
consuming  and  computer  intensive.  Thus,  further  analyses  in  this  area  of 
study  are  not  recommended  unless  a  very  specific  set  of  parameter 
conditions  needs  to  be  investigated.  However,  efforts  to  adapt 
continuation  software  tj  study  the  uncontrolled  equations  of  motion  would 
be  beneficial  tc  o  exploration  of  the  many  conservative  bifurcation 
points  found  »s  parameters  were  varied. 


Appendix  A.  Equations  of  Motion 


In  this  appendix,  the  equations  of  aotion  of  a  spinning  syMietric 
satellite  in  an  elliptical  orbit  about  a  spherical  attracting  body  are 
developed  via  Lagrange’s  nethod.  The  spin  axis  of  the  satellite  is  chosen 
to  be  in  the  plane  of  symaetry  of  the  satellite  and  perpendicular  to  the 
orbital  plane.  The  full  nonlinear  equations  of  motion  are  presented  along 
with  the  development  of  the  Lagrangian  and  the  body  angular  velocities. 

Kane  and  Barba  studied  the  stability  of  the  linearizred  equations  of 
motion  developed  herein  [8:402-405}.  An  alternative  development  for  the 
linearized  equations  of  motion  using  Euler’s  moment  equations  is  presented 
in  their  paper.  Captain  Dale  Shell  [13:2.23-2.25]  and  Captain  Jim  Cole 
[5:A-1-A-16]  presented  the  same  Newtonian  development  in  their  theses 
also. 

4.1  Attitude  Angular  Velocities 

The  following  develpoment  for  the  angular  velocites  parallels  that 
presented  by  Captain  Jin  Cole  [5:A-4-A-9].  However,  the  resulting  angular 
velocites  here  are  used  to  form  the  kinetic  energy  of  the  satellite  which 
I  will  call  T.  A  1-2-3  rotation  through  the  angles  0j,  Oj,  and  0j  measures 
the  orientation  of  the  satellite  with  respect  to  the  orbital  reference 
frame  A.  The  resulting  coordinate  frame  at  the  end  of  the  three 
consecutive  rotations  is  the  body  axis  frame  D.  Figure  A-1  shows  this 
graphically. 

Figure  A-2  shows  the  orientation  of  the  orbital  reference  frame  A 
with  repsect  to  an  inertial  reference  frame  I  which  has  its  origin  at  the 
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Figure  A.l  Coordinate  Reference  Fraaes 


center  of  the  attracting  body.  Axis  Aj  is  in  plane  of  the  orbit  and  points 
along  the  radius  vector  to  the  satellite  center  of  aass.  Axis  A^  lies  in 
the  orbital  plane  in  the  direction  of  the  aotion  of  the  satellite.  And 
finally,  Aj,  which  is  perpendicular  to  the  orbit  plane,  is  siaply  AjxAj  to 
fora  a  dextral,  orthogonal  coordinate  systea. 

Since  the  satellite  is  a  syaaetric  body,  we  can  siaplify  the 
developaent  of  the  expressions  for  the  angular  velocities  by  using  the 
nodal  coordinate  fraae  C.  The  nodal  fraae  is  achieved  after  the  second 
rotation.  Note  that  the  final  rotation  is  about  an  axis  of  syaaetry. 
Thus,  we  need  to  express  the  angular  velocity  of  the  body  fraae  D  attached 
to  the  satellite  with  respect  to  the  inertial  coordinate  fraae  I  written 


Figure  A. 2  Satellite  Reference  Fra^e  and  Orbit  Description 


in  nodal  coordinates  C.  This  is  sisply  the  sub  of  the  indivdual  rotation 
rates  to  go  froB  the  A  fraBe  to  the  D  fraBe. 

flO/xj  p+ac/*!  ^  ( A.  1 ) 

Now  let  us  look  at  each  individual  rotation  to  describe  its  rotation 
rate  and  coodinate  transforaBtion.  The  angular  velocity  of  the  A  fraae 
with  respect  to  the  inertial  fraae  is  the  rate  of  change  of  the  true 
anoBaly  for  the  orbit. 
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A 


(A. 2) 


The  transformation  matrix  for  this  rotation  is  not  needed.  The  other 
relative  rotation  rates  and  their  associated  coodinate  transformation 
matrices  are 


HB/a, 


B 


(A. 3) 


L 


BA 


1  0 
0  cosOi 
0  -aindi 


0 

sindj 

coa0i 


coa6,  0  -ainOj 
0  10 
ainO,  0  coaO, 


(A. 4) 


(A. 5) 


(A. 6) 


(A. 7) 


0080,  alnO,  0 

-alnOj  CO863  0 


0  0  1 


(A. 8) 


In  order  to  form  the  sum  of  these  components  of  the  angular 
velocity,  they  all  must  be  expressed  in  the  nodal  coordinate  frame  C.  To 
accomplish  this,  the  required  transformations  must  be  formed  such  that 
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c 


(A. 9) 


In  other  words,  we  aust  rewrite  equation  (A.l)  as 

flo/ii  ^.ac/flj  p  (A.  10) 

and  perfora  the  aatrix  aultiplication  to  transfora  all  the  angular 
velocities  to  the  C  fraae.  The  transforaation  aatrices  needed  in  equation 
(A. 10)  are 


cosdj  0  -flinOj 
0  10 
alnS,  0  cosO, 


(A. 11) 


cosO, 

0 

slnO, 


sln6]^ain03 

CO801 

-8in0iCO80j 


-co80^8in02 

ain0]^ 

0080^0083 


(A. 12) 


is  not  needed  to  fora  the  angular  velocity,  but  it  will  be  needed  later 
for  the  derivation  of  the  potential  energy.  So,  it  is  given  as 


h  Iz  h 

®a 

k  ^ 


(A. 13) 


where 


ij^-CO803CO8O, 

(A. 14) 

i,  -  8  InO^a  ln03CO80,  co80j^8  in03 

(A. 15) 

i3--co8O38ln03CO80,<t-8in0]^8ln0, 

(A. 16) 
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/n^--coB028in03 

(A. 17) 

in2--sin02sln02ain03 -^0080300803 

(A. 18) 

in3-oo8O38in02sln03^8in03OO803 

(A. 19) 

7l3-8ln02 

(A. 20) 

ll2“-8in0iOO80j 

(A. 21) 

713-0080300802 

(A. 22) 

Substituting  the  transforaations  into  equation  (A. 10)  and  performing 
the  matrix  algebra  results  in  the  expessions  for  the  angular  velocity  of 
the  satellite  expessed  in  nodal  coordinates. 

[  djCosOj-^cosOiSlndj 

e^+^aine^  (A.  23) 

P3+0lSin02+^CO8eiCO862  ^ 

\.2  Kinetic  Energy  of  the  Satellite 

The  expression  for  the  kinetic  enregy,  T,  is  needed  to  fora  the 
Lagangian  which  is  used  to  derive  the  equations  of  motion.  A  good 
presentation  of  this  derivation  is  in  Wiesel  [17:103-104]  and  is  presented 
here  for  completeness.  The  general  expression  for  the  T  is 

T'*  J (A. 24) 

Figure  A-2  shows  the  position  vector  R  to  the  center  of  mass  of  the 
satellite  and  the  vector  p  to  a  mass  element  da  in  the  body.  The  velocity 
of  the  mass  element  da  is 
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(A. 25) 


where 


dt  dt 


(A. 26) 


Substituting  these  expessions  into  equation  (A. 24)  results  in  the 
following  integral  over  the  body. 


(A. 27) 


But,  since  the  sass  eleaent  da  is  Measured  from  the  center  of  aass, 


(A. 29) 


So,  equation  (A. 28)  reduces  to 

r=-|i»-^-^‘-i/(8‘’"xii»  J  •  (8‘’"x{?)  J  dm  (A.30) 


where 


(A. 31) 


Substituting  equation  (A. 31)  into  equation  (A.30),  we  get 
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(A. 32) 


+ (i)?(x3  +X3 )- 2 WjWjXjXj  -2«i«3XjX3 - 2 }  <in 
Froa  rigid  body  dynamics,  the  moment  of  inertia  matrix  can  be  written  as 

j’(xa+x|)dm  -jx^x^dm 

J»  -jx^x^dm  j{xl+xf)dm  -jx^dm  (A. 33) 

-JxiX3<im  -Jj^X3<±n  j{xl*xj)da 

Let  us  use  a  principal  body  axes  frame  to  simplify  the  moment  of  inertia 
matrix  and  let  the  principal  moments  of  inertia  be  A,  B  ,  and  C.  For 
further  reading  on  this  subject  see  reference  [17:102].  For  a  symmetric 
body,  B=A,  thus  the  moment  of  inertia  matrix  becomes 

A  0  0]  [A  0  0 

J-  0  B  0  ■  0  A  0  (A. 34) 

p  0  Cj  [o  0  C 

Substituting  the  moment  of  inertia  terms  into  equation  (A. 32)  reduces  the 
expression  for  the  kinetic  energy  for  the  satellite  to 

This  is  simply  the  kinetic  energy  of  the  translational  notion  of  the 
center  of  mass  plus  the  kinetic  energy  of  the  rotational  motion  about  the 
center  of  mass.  The  velocity  of  the  center  of  mass  is 

After  performing  the  cross  product,  we  get 
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Since  the  satellite  is  a  synaetric  body,  we  can  use  the  angular  velocity 
written  in  the  nodal  coordinates  to  siaplify  the  calculation  of  the 
kinetic  energy.  Substituting  the  angular  velocity  components  of  equation 
(A. 23)  and  the  velocity  of  the  center  of  mass  above,  we  get  an  expression 
for  the  kinetic  energy  that  looks  like 


T=  -iA(diCOs02-^cos6j^8ln0j)* 

2k  2 

*  A(0j + a  in0  C(03  +0^  s  in02 + ♦  cosO^co  s02)* 

2  2 


(A. 38) 


A, 3  Potential  Energy  of  the  Satellite 

The  potential  energy  of  the  satellite,  V,  is  the  gravitational 
potential.  This  is  derived  in  Meirovitch  [9:430:438]  and  presented  here 
in  further  detail.  In  order  to  derive  the  potential,  we  Bust  locate  a 
aass  eleaent  da  in  the  satellite  body  with  respect  to  the  inertial  fraae 
at  the  origin  of  the  attracting  body.  Figure  A-2  shows  the  vector  B  from 
the  center  of  the  attracting  body  to  the  center  of  aass  of  the  satellite, 
the  vector  p  froa  the  satellite  center  of  aass  to  the  mass  element  da,  and 
the  vector  r  such  that 


The  potential  can  be  written  as 


{A. 39) 


(A. 40) 


where  G  is  the  universal  gravitational  constmt  and  M  is  the  aass  of  the 


attracting  body.  We  can  write  the  magnitude  of  the  vector  r  as 


(A. 41) 

We  can  factor  out  an  R  so  that 

where  e^  is  the  unit  vector  in  the  direction  of  the  R  vector.  Taking  the 
inverse  of  this  expession  results  in 


which  allows  us  to  readily  expand  it  via  the  binomial  expansion  keeping 
only  second  order  terms. 


(1  +  e  )•*  1  1 ) 


2t 


+  ••• 


The  result  after  the  expansion  is 


(A. 44) 


r 


(A. 45) 


Substituting  the  expansion  into  equation  (A. 40)  for  the  potential 
calculation,  we  get 


Again,  since  p  is  measured  from  the  mass  center, 

S,  f 

— -•  I  ddn»0 


and  equation  (A. 46)  reduces  to 


(A. 46) 


(A. 47) 
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(A. 48) 


The  unit  vector  in  the  radial  direction  can  be  expressed  in  the  body  fraae 
by  the  following: 


(  a.  49 ) 

Let  us  define  the  vector  p  in  the  body  fraae  D  to  be 

(A.  50) 


Substituting  equations  (A. 49)  and  (A. 50)  into  the  reaaining  integral  of 
equaution  (A. 48)  leaves  the  potential  as 


*(3iii  - ijxj  *62 ^m^x^x^*6m^n^x^-^*61  }  dm 


(A. 51) 


With  soae  aanipulation  of  the  integrals  in  the  eoaent  of  inertia  Matrix  of 
equation  (A. 33),  the  following  integrals  can  be  shown  to  be 

fxidaf-^(A*B*Q  (A.  52) 

f  Xa<im»~(A+C-fl)  (A. 53) 

f^XjdOf^iA+B-  O  { A .  54  ) 

We  again  assuae  a  set  of  principal  axes  to  null  the  integrals  of  the 
product  tens.  Making  the  appropriate  substitutions  into  equation  (A. 51) 
and  perforaing  the  algebra,  the  potential  reduces  to 
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v^- -  -  l^|ii^(B+C-A)+ffli^(A+C’-B)+/ii*(;i+B-C)- -|•(-A+B+Q|  < A.  55 ) 

Substituting  in  for  1.,  ij ,  and  nj  and  Baking  B=A  in  the  aoaent  of  inertia 
aatrix,  the  potential  becoaes 

.^^Ctoa^2+2Aflin*8j-Cbin*02--j(2A+Q|  (A.  56) 

But,  using  the  tr igonoaetric  identity 

aBin^2  +  0CO8^6j-C  (A.  57) 

we  can  show  that  the  potential  siaply  becoaes 

V-  -^(A-09in>e,  -  -j(A-o]  ( A.  58 ) 


A. 4  Lagranges  Equations  of  Motion 

Lagrange’s  equations  of  aotion  for  holonoaic,  conservative, 
dynaaical  systeas  are  foraed  by  using  the  equation 


d  dL  dL 
dt  d<t^  dq^ 


k*l,2, '”,a 


(A. 59) 


An  excellent  reference  on  this  aaterial  is  Meirovitch’s  text  on  analytical 
dynaaics.  [9]  The  q^’s  are  the  generalized  coordinates  for  the  systea  and 
L  is  teraed  the  Lagrangian  which  is  siaply 


L»T-V 


(A. 60) 


where  T  is  the  kinetic  energy  and  V  is  the  potential  energy.  These  are 
given  respectively  by  equations  (A. 38)  and  (A. 58)  froa  the  previous 
derivations.  So,  substituting  T  and  V  gives  the  Lagrangian 
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A  A  «  VA»01/ 

♦  -  j  C(6  3  •*•02  8  ln02  *  ^  cobQ^cobQ^Y  *  ♦  -  j  |(A-  0^8  in*0j  -  -i 

The  nice  thing  about  analytical  dynamics  is  that  constants  of  the  notion 
are  identified  rather  easily.  In  this  case,  one  of  the  constants  of 
notion  is  identified  by  the  coordinate  Sj  not  appearing  in  the  Lagrangian, 
thus  naking  it  an  ignorable  coordinate.  Further,  since 


dL 

00, 


(A. 62) 


then  it  follows  that  the  conjugate  generalized  angular  nonentun. 


dL 


(A. 63) 


is  constant  by  rearrangement  of  equation  (A. 59)  to  show  that 


dL  ^  dL 


dt  a0,  "557 


(A. 64) 


Note  that  the  coordinate  v  does  not  appear  in  the  Lagrangian  either. 
However,  we  will  uses  this  constant  of  the  notion  later  to  get  one  of  the 
equations  of  notion.  Also,  I  would  like  to  nention  that  the  Lagrangian 
does  not  contain  tine  explicitly  which  leads  to  a  constant  of  the  notion 
known  as  Jacobi’s  integral.  However,  we  will  not  need  Jacobi’s  integral 
to  forn  our  equations  of  notion  so  it  will  not  be  presented  here.  On  the 
other  hand,  it  is  worthwhile  to  nention  that  the  Jacobi’s  integral  is  a 
statenent  of  the  law  of  conservation  of  energy. 

Since  the  Lagrangian  for  the  satellite  dynanics  has  two  ignorable 
coordinates  and  does  not  contain  tine  explicitly,  we  can  use  Routh’s 
method  for  the  ignoration  of  coordinates.  However,  we  will  only  use  this 
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technique  to  reiaove  the  constant  of  motion  for  the  coordinate  6,.  For  a 
detailed  explanation  of  this  method,  you  are  again  referred  to 
Meirovitch’s  text.  [9:85-88]  The  Routhian  function  is 

R=L-P^t^  (A.  65) 

where 

+^^8  in0j coaBjCosOj)  { A .  66 ) 

After  making  the  appropriate  substitutions  for  L  and  some  work,  the 
Routhian  function  becomes 

-|■A(6lCO80J-^cos^^8ln0J)* 

M  A 

♦  •jA(0j+^Sin0i)^+Pj^(0l8ln0j+^CO86^COB02)  (A.  67) 

-  4)] 

It  can  be  shown  that  the  equations  of  motion  for  the  remaining  coordinates 
are  formed  by  using  equation  (A. 59)  and  replacing  L  with  R.  The 
solutions  to  these  equations  allow  us  to  solve  for  the  ignorable 
coordinate,  0j ,  using  the  constant  of  motion  given  by  equation  (A. 63). 

We  are  now  ready  to  form  the  equations  of  motion  using  R,  0j,  and 
02  as  our  generalized  coordinates.  Let’s  start  by  deriving  the  equation 
for  the  coordinate  ?.  The  partial  derivatives  and  the  derivative  with 
respect  to  time  are 

(A. 68) 

-4rP^^mR^^*2aRjt^  (A.  69) 

oft 
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(A. 70) 


I 


iE. 

dv 


»0 


Note  that  we  assumed  that  the  attitude  motion  of  the  satellite  does  not 
effect  the  orbital  motion,  thus  the  approximation  for  given  by  equation 
(A. 68).  The  coordinate  v  does  not  appear  in  the  Lagrangian,  so  is  a 
constant  of  motion  also  which  merely  states  that  the  angular  momentum  of 
the  orbit  is  conserved.  We  now  write  the  specific  angular  momentum,  H,  as 
a  function  of  the  two  body  orbital  elements 


(A. 71) 


where  the  mean  motion  of  the  satellite  is 


/J- 


2H 

T 


(A. 72) 


and  the  gravitational  parameter  is  approximated  by 

(A. 73) 

Any  fundamental  text  in  astrodynamics  such  as  Bate,  Mueller,  and  White 
shows  this  in  further  detail  [2).  So,  our  equation  for  the  coordinate  v 
is 

(A. 74) 

The  derivatives  needed  to  for  the  equation  of  motion  for  the 
coordinate  R  are 

(A. 75) 

die 


A-15 


(A. 76) 


_d_ 

dt 


dR  GMa 


(A. 77) 


Again,  we  assume  that  the  attitude  motion  of  the  satellite  does  not  effect 
the  orbit  by  the  approximation  in  equation  (A. 77).  So,  dividing  by  m  and 
using  the  relation  in  equation  (A. 73),  we  get 

(A. 78) 

R^ 

or 


(A.  79) 

R^ 

Now  we  substitute  equation  (A. 74)  into  the  above  to  get  the  equation  for 
B  in  terms  of  the  orbital  elements. 


R^  J?* 


(A. 80) 


It  is  to  our  advantage  to  nondimensionalize  by  introducing  the  rew 
variables 


T»nt 


(A. 81) 


(A. 82) 


with 


(A.83) 

aZ 

In  order  to  nondimensionalize  equations  (A. 74)  and  (A. 79),  we  need  to  find 
all  derivatives  with  respect  to  t,  which  will  be  represented  by  primes. 
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So,  applying  the  chain  rule,  we  get 


and 


using  our  new  variable  definitions.  If  we  divide  equation  (A. 80) 
we  get 


(A 


{A. 85) 
by  n  a, 


4. -gl  =  0  (A. 86) 

n^a 

Then  by  using  equations  (A. 84)  and  (A. 85),  we  get  the  transformed  equation 
for  the  R  coordinate  to  be 


'  *2 


(A. 87) 


C-  V- 

To  nondimensional ize  equation  (A. 74),  we  use  the  chain  rule  again  to  get 


^  cfv  dt  _  a*r(X-0^)  ^  1 
dt  dt  dx  j?2  n 


which  becomes 


(A. 88) 


(A. 89) 


We  now  take  the  derivative  with  respect  to  r  to  get 


-2(i-»»)'^  r, 

C’ 


(A. 90) 


which  we  will  need  later. 

Now  that  we  have  the  nondimensionalized  orbital  equations,  we  must 
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make  the  appropriate  coordinate  transformation  on  the  initial  conditions. 
We  arbitrarily  select  the  satellite  to  be  at  perigee  at  t=T=0.  So,  in  the 
old  coordinates, 

/a(0)»a(l-e)  (A.  91) 

and 

(A.  92) 

are  the  initial  conditions.  In  the  new  nondimensional  coordinates,  we  get 

C(0)»l-O  (A.  93) 

and 

(A.  94) 

for  our  initial  conditions. 

Next,  we  turn  our  attention  to  the  attitude  equations  of  motion. 
The  required  partial  and  tine  derivatives  for  the  9j  coordinate  are 

Pe^=-^*-A(diCO803-^cos6j^Bin82)co862+P«  siiiOj  (A.95 ) 

Pj^^0jCC8^6j -2ii6263CO362CO80j8in6j 

-AV02CO802(cO8*6j-8in*02)+At0i8in0i8in0jCO802  ( A.  96 ) 

-il^CO80j^CO8028in02^P^d2<^^®2 


%I(02CO802 -♦cosBiBlnBjjffr  BlnBjSlnBi) 
+A(02  ♦♦  8  in0i)(tco86i)  -  8  in02CO862 


(A. 97) 


After  substituting  the  derivatives  into  equation  (A. 59),  Lagrange’s 
equation,  we  get 


A- 18 


-<lcos63^tan02+-^d28ec62-^di8in0itan62 

+ <i^cos0iS  in0i  t  an^02  -62^  CO  a0i(l  ♦  tan*02) 

p 

-^*cos0isin0i(l  +  tan^02)+  — sinOiSec02 


(A. 98) 


V/e  see  that  six  teras  conveniently  cancel  out  and  two  terms  combine  to 
simplify  to 


Finally, 


O«6^-20^02tan02-2V02CO802-S>co80item02 


-V*co80i8in0i+— |i8ec02(d2+^8in0i) 


(A. 99) 


02 ’■20^02  tan62  +  2^03COS02+^COS0^t2U102 

p 

♦♦*co802Bin0i-(jr+l)-^8ec02(02+'>8in02) 

is  the  equation  of  motion  for  the  6|  coordinate  where 

C-A 


(A. 100) 


(A. 101) 


and 


C 


(A. 102) 


is  the  constant  of  motion  which  is  the  spin  rate  about  the  symmetry  axis. 
The  required  derivatives  for  the  Oj  coordinate  are 


i’a,*-^»A(02+'»8in02) 


(A. 103) 


■—  »A02  +^^0200802  ♦AiJ  8  in02 


(A. 104) 
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(A. 105) 


-  >sA(diCos6j  -♦cos0iain62)(-di8in6j-^co88iCO802) 

I 

+P^(0ico302-Vcos0i8in02)+^^(A-C)8in02cos02 


Multiplying  out  the  teras  reduces  equation  (A. 105)  to 

A0?8in02CO802-A0iVcos02CO8*02+A^03^co80i8in*03 

oOj 

+Ai*co8^j^8in02Coa02+P#j(ftiCO802-^co80i8in02) 

+  -^(A-  C)8in02CO802 
R 


(A. 106) 


Substituting  the  derivatives  into  Lagrange’s  equation  one  last  time 
results  in 

-0lVCO80jL8in*02-^^CO8*0l8ln02CO^2  (A  107) 

— ^(0iCO802-’>co80iainO2)+-^|-~^j8in0jCO802 

for  the  Oj  coordinate.  Rearranging  the  above  equation  and  substituting  in 
for  K  reduces  the  equation  of  aotion  to 


(0^00802 -'>CO8028in02) 

-3  Aain0aCO802-f  8in0, 

^3  a  a 

-0j8in02COa02+^*CO8*028in02CO803 
-♦0jCO8Oi  ♦V02CO80i(8in*0a -cos^Oj) 


(A. 108) 


And  finally,  using  the  trigonoaetric  identity 


8  in^2  ■  c®®  *®2  “  ^  ■  2co®  *®2 


(A. 109) 


the  equation  of  aotion  for  the  02  coordinate  siaplifies  to 
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fl  j  ■(#:+ 1 ) -^  (6 1COS62  -  ^  coadjS  in6  j)  -  8  i  n6i 


(A, 110) 


Now  we  nondimensionalize  equations  (A. 100)  and  (A. 110)  using  the  new 
nondinensional  coordinates  given  by  equations  (A. 81)  and  (A. 82).  Using 
the  chain  rule,  it  is  easy  to  see  that 


c®!  .  c»i  dr 

dt  "  dx  dt  dx  ~  ^ 


(A. Ill) 


df<=®i)_  df  dr  ,__a 

dt(  dt  j  dri,  dr  j  dt  ^ 


(A. 112) 


are  obvious.  Using  these  results  with  equation  (A. 100)  lakes  the  0j 
coordinate  equation  of  aotion  transfon  to 

&(n  *  -  2  iJ  *0^0^  tan0j  ♦2i5  ’  v'0aCOB02 +n  *  v"co80i  tanO, 

«  •Piu//  .  .  (A. 113) 

♦jj*v'’cos028iii0j-(Jr+l)n*— ^(0i+v'8in0i)8ecO3 


where 


*  c(0'3  ln02  ♦  v^cobOjCobOj) = C(8+ 1) 


(A. 114) 


after  the  transforiation.  Now,  dividing  both  sides  by  n‘  gives  the 
transformed  result  for  the  9j  coordinate  as 

0$^- 2010^  tanOj  ♦  2  i'02CO802  ♦  v"co80i  tanOj  +  v^co80i  8  in0i 

P#  ,  .  ,  ,  (A. 115) 

-(K+ 1)  — f  {0^ + v'8  in02)aec03 


Next,  using  the  relations 
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(A. 116) 


and 


1  _ 

C"*  J?3 


(A. 117) 


we  reBove  the  R  froa  the  62  equation  of  aotion.  Then,  in  order  to 
transform  the  0,  equation  of  motion  to  the  nondimensional  coordinates,  we 
follow  the  exact  steps  as  we  did  for  the  0j  equation  of  motion.  The  final 
result  for  the  0,  coordinate  is 


®2-(^+l)-^(®lco8e2-v'coseiSineJ  (A.  118) 

-|0i^-v'^coa*0j-*--|^jaln02CO882-v"8in0i-2v^0icos6jCO8^2 


Equilibrium  conditions  of  the  attitude  equations  can  now  be 
determined.  The  equilibrium  condition  arises  when 


o'-0'3-e''-9"-o 


The  values  of 


(A. 119) 


01-03-0  (A. 120) 

malce  these  conditions  true.  In  this  case, 

03-(a^l)T-v-(a+i)ijt'W  (A.  121) 

This  equilibrium  condition  makes  the  satellite  spin  axis  perpendicular  to 
the  orbit  plane.  In  this  thesis,  this  is  the  equilibrium  the  satellite 
control  system  is  designed  about. 
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